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Abstract: We discuss several connections between discrete and continu- 
ous random trees. In the discrete setting, we focus on Galton- Watson trees 
under various conditionings. In particular, we present a simple approach 
to Aldous' theorem giving the convergence in distribution of the contour 
process of conditioned Galton- Watson trees towards the normalized Brow- 
nian excursion. We also briefly discuss applications to combinatorial trees. 
In the continuous setting, we use the formalism of real trees, which yields 
an elegant formulation of the convergence of rescaled discrete trees towards 
continuous objects. We explain the coding of real trees by functions, which 
is a continuous version of the well-known coding of discrete trees by Dyck 
paths. We pay special attention to random real trees coded by Brown- 
ian excursions, and in a particular we provide a simple derivation of the 
marginal distributions of the CRT. The last section is an introduction to 
the theory of the Brownian snake, which combines the genealogical struc- 
ture of random real trees with independent spatial motions. We introduce 
exit measures for the Brownian snake and we present some applications to 
a class of semilinear partial differential equations. 
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Introduction 

The main purposes of these notes are to present some recent work about con- 
tinuous genealogical structures and to point at some of their applications. The 
interest for these continuous branching structures first arose from their connec- 
tions with the measure-valued branching processes called superprocesses, which 
have been studied extensively since the end of the eighties. Independently of 
the theory of superprocesses, Aldous PP,[2] discussed scaling limits of various 
classes of discrete trees conditioned to be large. In the case of a Galton- Watson 
tree with a finite variance critical offspring distribution and conditioned to have 
a large number of vertices, he proved that the scaling limit is a continuous ran- 
dom tree called the Brownian CRT. Moreover, this limiting continuous object 
can be coded by a normalized Brownian excursion, a fact that is reminiscent 
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of the Brownian snake construction of superprocesses |25| . In the recent years, 
these ideas have been extended to much more general continuous trees: See in 
particular ^5] for a discussion of Levy trees, which are the possible scaling limits 
of critical or subcritical Galton- Watson trees. 

Section below is concerned with scaling limits of discrete trees. In fact, we 
do not really discuss limits of trees, but consider rather certain coding functions 
of these trees, namely the height function and the contour function (see Fig^ 
below). Our main result ( Theorem 1 1.8fl states that the rescaled height function 
associated with a forest of independent (critical, finite variance) Galton- Watson 
trees converges in distribution towards reflected Brownian motion on the pos- 
itive half-line. From this, one can derive a variety of analogous limit theorems 
for a single Galton- Watson tree conditioned to be large. The derivation is quite 
simple if the conditioning is "non-degenerate" : For instance if the tree is condi- 
tioned to have height greater than n (resp. total progeny greater than n), the 
scaling limit of the height function will be a Brownian excursion conditioned 
to have height greater than 1 (resp. duration greater than 1). For degenerate 
conditionings, things become a little more complicated: The case of a Galton- 
Watson tree conditioned to have exactly n vertices, corresponding to Aldous 
theorem |2j, is treated under an exponential moment assumption using an idea 
of Marckert and Mokkadem |3J . We briefly discuss some applications to various 
classes of "combinatorial" trees. 

Although the limit theorems of Section^give a lot of useful information about 
asymptotics of discrete random trees, it is a bit unsatisfactory that they only 
discuss continuous limits for the coding functions of trees and not for the trees 
themselves. The formalism of real trees, which is briefly presented in Section 
provides an elegant way of restating the limit theorems of Section ^ in terms 
of convergence of trees. The use of real trees for probabilistic purposes seems 
to be quite recent: See in particular ^2] and In Section we first discuss 
the coding of a (compact) real tree by a continuous function. This is of course 
a continuous analogue of the correspondence between a discrete tree and its 
contour function. This coding makes it possible to get a simple and efficient 
construction of the CRT and related random trees as trees coded by various 
kinds of Brownian excursions. As an application, we use some tools of Brownian 
excursion theory to derive the finite-dimensional marginals of the CRT, which 
had been computed by Aldous with a very different method. 

Section |31 gives an introduction to the path- valued process called the Brow- 
nian snake and its connections with certain semilinear partial differential equa- 
tions. The Brownian snakes combines the genealogical structure of the random 
real trees studied in Section|21with spatial motions governed by a general Markov 
process £. Informally, each Brownian snake path corresponds to the spatial posi- 
tions along the ancestral line of a vertex in the tree. The precise definition of the 
Brownian snake is thus motivated by the coding of real trees that is discussed in 
Section[21 In view of applications to PDE, we introduce the exit measure from a 
domain D, which is in a sense uniformly spread over the set of exit points of the 
Brownian snake paths from D. We then derive the key integral equation (Theo- 
rem f° r the Laplace functional of the exit measure. In the particular case 
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when the underlying spatial motion £ is <i-dimensional Brownian motion, this 
quickly leads to the connection between the Brownian snake and the semilinear 
PDE Au = u 2 . Up to some point, this connection can be viewed as a reformula- 
tion of (a special case of) Dynkin's work J2| about superprocesses. Although we 
content ourselves with a simple application to solutions with boundary blow-up, 
the results of Section |21 provide most of the background that is necessary to un- 
derstand the recent deep applications of the Brownian snake to the classification 
and probabilistic representation of solutions of Au — u 2 in a domain |33| . 

The concepts and results that are presented here have many other recent 
applications. Let us list a few of these. Random continuous trees can be used to 
model the genealogy of self-similar fragmentations [20] ■ The Brownian snake has 
turned out to be a powerful tool in the study of super-Brownian motion: See the 
monograph |27| and the references therein. The random measure known as ISE 
(see [3] and Section [3] below), which is easily obtained from the Brownian snake 
driven by a normalized Brownian excursion, has appeared in asymptotics for var- 
ious models of statistical mechanics [J] , [23 > [2H] • There are similar limit theorems 
for interacting particle systems in X d : See |S] for the voter model and coalescing 
random walks and the recent paper (23 f° r the contact process. Another promis- 
ing area of application is the asymptotic study of planar maps. Using a bijection 
between quadrangulations and well-labelled discrete trees, Chassaing and Scha- 
effer were able to derive precise asymptotics for random quadrangulations in 
terms of a one-dimensional Brownian snake conditioned to remain positive (see 
PH) for the definition of this object). The Chassaing-Schacffer asymptotics have 
been extended to more general planar maps by Marckert and Miermont \62\. 
In fact one expects the existence of a universal continuous limit of planar maps 
that should be described by random real trees of the type considered here. 

1. Prom Discrete to Continuous Trees 

In this section, we first explain how discrete random trees can be coded by dis- 
crete paths called the height function and the contour function of the tree. We 
then prove that the rescaled height function associated with a forest of indepen- 
dent Galton- Watson trees converges in distribution towards reflecting Brownian 
motion on the positive half-line. This has several interesting consequences for 
the asymptotic behavior of various functionals of Galton- Watson forests or trees. 
We also discuss analogous results for a single Galton- Watson tree conditioned 
to be large. In particular we recover a famous theorem of Aldous showing that 
the suitably rescaled contour function of a Galton- Watson tree conditioned to 
have n vertices converges in distribution towards the normalized Brownian ex- 
cursion as 7i — > oo. Consequences for various classes of "combinatorial trees" are 
outlined in subsection 1 1.51 
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1.1. Discrete trees 

We will be interested in (finite) rooted ordered trees, which are also called plane 
trees in combinatorics (see e.g. [H]). We first introduce the set of labels 

oo 

w=|Jr 

n=0 

where N = {1,2, . . .} and by convention N° = {0}. An element of U is thus 
a sequence u — (it 1 , . . . , u n ) of elements of N, and we set |u| = n, so that |u| 
represents the "generation" of u. If it = (it 1 , . . . u m ) and v = (v 1 , . . . , v n ) belong 
to U, we write uv = (it 1 , . . . u m , v 1 , . . . , v n ) for the concatenation of u and v. In 
particular u0 — 0u — u. 

The mapping it : U\{0} — > U is defined by ^(u 1 . . . u n ) — u 1 . . . u n ~ l (w(u) 
is the "father" of u). 

A (finite) rooted ordered tree t is a finite subset of U such that: 

(i) e t. 

(ii) u S t\{0} ^> n(u) G t. 

(iii) For every wet, there exists an integer k u (t) > such that, for every 
j G N, uj G t if and only if 1 < j < k u (t) 

The number fc u (t) is interpreted as the "number of children" of u in t. 

We denote by A the set of all rooted ordered trees. In what follows, we see 
each vertex of the tree t as an individual of a population whose t is the family 
tree. The cardinality #(t) of t is the total progeny. 




Tree t Contour function C 3 Height function ht(n) 

Figure 1 

We will now explain how trees can be coded by discrete functions. We first 
introduce the (discrete) height function associated with a tree t. Let us denote 
by Mo = 0,iti, U2, . . . , it#(t)-i the elements of t listed in lexicographical order. 
The height function (ht(n); < n < #(t)) is defined by 

ftt(n) = |«n| ) 0<n<#(t). 
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The height function is thus the sequence of the generations of the individuals 
of t, when these individuals are listed in the lexicographical order (see Fig^for 
an example). It is easy to check that ht characterizes the tree t. 

The contour function (or Dyck path in the terminology of |¥T|) gives an- 
other way of characterizing the tree, which is easier to visualize on a picture 
(see Fignj. Suppose that the tree is embedded in the half- plane in such a way 
that edges have length one. Informally, we imagine the motion of a particle that 
starts at time t = from the root of the tree and then explores the tree from 
the left to the right, moving continuously along the edges at unit speed (in the 
way explained by the arrows of Fig^), until all edges have been explored and 
the particle has come back to the root. Since it is clear that each edge will 
be crossed twice in this evolution, the total time needed to explore the tree is 
C(t) := 2(=/^(t) — 1). The value C s of the contour function at time s £ [0, C(t)] 
is the distance (on the tree) between the position of the particle at time s and 
the root. By convention C s = if s > C(t). Fig^ explains the construction of 
the contour function better than a formal definition. 

We will introduce still another way of coding the tree. We denote by S the 
set of all finite sequences of nonnegative integers mi, ... , m p (with p > 1) such 
that 

• mi + m 2 + • • • + rrii > i , Vi € {1, . . . ,p — 1}; 

• mi + m,2 -\ 1- m p = p — 1. 

Recall that u a — 0, u%, U2, ■ ■ ■ , w#(t)-i are the elements of t listed in lexico- 
graphical order. 

Proposition 1.1 The mapping 



defines a bijection from A onto S. 

Proof. We note that if #(t) = p, the sum k UQ (t)+fc Ul (t)+- • •+fc u#(t) _ 1 (t) counts 
the total number of children of all individuals in the tree and is thus equal to p— 1 
(because is not counted !). Furthermore, if i g {0, 1, . . . ,p— 2}, k UQ +- ■ -+k Ui is 
the number of children of uo, ... ,Ui and thus greater than or equal to i, because 
Ui,. .. ,Ui are counted among these children (in the lexicographical order, an 
individual is visited before his children) . There is even a strict inequality because 
the father of Ui+i belongs to {«o, . . . , Ui}. It follows that $ maps A into S. We 
leave the rest of the proof to the reader. □ 

Let t 6 A and p — #(t). Rather than the sequence (mi, . . . , m p ) — $(t), we 
will often consider the finite sequence of integers 



$ : t 



{k Uo (t): k Ul (t), . . . 



(t)) 



n 




< n < p 



which satisfies the following properties 
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• so = and x p = — 1. 

• x n > for every < n < p — 1 . 

• Xi — Xi-i > — 1 for every 1 < i <p. 

Such a sequence is called a Lukasiewicz path. Obviously the mapping $ of the 
proposition induces a bijection between trees (in A) and Lukasiewicz paths. 

We now observe that there is a simple relation between the height function 
of a tree and its Lukasiewicz path. 

Proposition 1.2 The height function ht of a treet is related to the Lukasiewicz 
path of t by the formula 

ht(n) = #{i S {0, 1, . ..,n - 1} : xj = inf xg}, 

j<£<n 

for every n € {0, 1, . . . ,#(t) - 1}. 
Proof. Obviously, 

ht(n) = #{j 6 {0, 1, . . . ,n - 1} : Uj -< u n } 

where -< stands for the genealogical order on the tree (it -< v if v is a descendant 
of u). Thus it is enough to prove that for j £ {0, 1, . . . , n — 1}, Uj <u n holds iff 

Xj = inf Xe. 

j<£<n 

To this end, it suffices to verify that 

inf{fc > j : Xk < Xj} 

is equal either to #(t), in the case when all Uk with k > j are descendants of 
Uj, or to the first index k > j such that is not a descendant of Uj. 
However, writing 

k 

Xk - Xj = ^ ( m i _ !) 
i=j+l 

and using the same arguments as in the proof of Proposition ll.il (to prove that 
$ takes values in S), we see that for every I > j such that ug is a descendant 
of Uj, we have xg — Xj > 0, whereas on the other hand Xk — Xj = —1 if k is the 
first I > j such that ug is not a descendant of j (or fc = p if there are no such 
£). This completes the proof. □ 



1.2. Galton-Watson trees 

Let fi be a critical or subcritical offspring distribution. This means that \i is a 
probability measure on Z + such that 

oo 

£m*)<i- 

fc=0 
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We exclude the trivial case where /i(l) = 1. 

We will make use of the following explicit construction of Galton- Watson 
trees: Let (K u , u G hi) be a collection of independent random variables with law 
/x, indexed by the label set hi. Denote by 9 the random subset of hi defined by 

9 = {u = u 1 . . . u n e hi : v? < K u x_ M 3-i for every 1 < j < n}. 

Proposition 1.3 9 is a.s. a tree. Moreover, if 

Z n = #{ue9: \u\ = »}, 

(Z n ,n > 0) is a Galton-Watson process with offspring distribution fi and initial 
value Zq = 1. 

Remark. Clearly k u (9) = K u for every u G 9. 

The tree 9, or any random tree with the same distribution, will be called a 
Galton-Watson tree with offspring distribution /i, or in short a /i-Galton- Watson 
tree. We also write n M for the distribution of 9 on the space A. 

We leave the easy proof of the proposition to the reader. The finiteness of 
the tree 9 comes from the fact that the Galton-Watson process with offspring 
distribution /x becomes extinct a.s., so that Z n = for n large. 

If t is a tree and 1 < j < fc (t), we write Tjt for the tree t shifted at j: 

Tjt = {ueU : ju e t}. 

Note that Tjt is a tree. 

Then n p may be characterized by the following two properties (see e.g. [Hlj 
for more general statements): 

(i) U ll (k = j) = j e Z+. 

(ii) For every j > 1 with fi(j) > 0, the shifted trees Tit, . . . , Tjt are independent 

under the conditional probability H^(dt \ k — j) and their conditional 
distribution is 11^. 

Property (ii) is often called the branching property of the Galton-Watson 
tree. 

We now give an explicit formula for 11^. 
Proposition 1.4 For every t G A, 

n,(t) = H^(t)). 

net 

Proof. We can easily check that 

{9 = t}= f]{K u = k u (t)}, 

net 
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so that 

n M (t) =P(6 = t) = Y[ P(K U = fc„(t)) = J] /i(fc„(t)). 

□ 

Recall from Proposition II . II the definition of the mapping <I>. 
Proposition 1.5 Let 9 be a /i-Galton- Watson tree. Then 

t>(6) ^ (M U M 2 ,...,M T ), 

where the random variables M\, M 2 , ■ ■ ■ are independent with distribution \x, and 

T = inf{n > 1 : M x H V M „ < n}. 

Remark. The fact that T < oo a.s. is indeed a consequence of our approach, 
but is also easy to prove directly by a martingale argument. 

Proof. We may assume that 9 is given by the preceding explicit construction. 
Write Uo = 0, U\, . . . , £/#(g)_i for the elements of 9 listed in lexicographical 
order, in such a way that 

$(6) = (Ku ,K Ul ,...,K u#w _ 1 ). 

We already know that K Uo + ■ ■ - + K Un > n+ 1 for every n £ {0, 1, ... , #(0) — 2}, 
and K Uo + ■ ■ ■ + K u#m _ t = #(0) - 1. 

It will be convenient to also define U p for p > #(0), for instance by setting 

Up = E/#(0)_il ... 1 

where in the right-hand side we have added p — #(0) + 1 labels 1. Then the 
proof of the proposition reduces to checking that, for every p > 0, Kjj , . . . , Kjj 
are independent with distribution /i. The point is that the labels Uj are random 
(they depend on the collection (K u , u € U)) and so we cannot just use the fact 
that the variables K u , u £ U are i.i.d. with distribution /x. We argue by induction 
on p. For p = or p = 1, the result is obvious since Uq = and U± = 1 are 
deterministic. 

Fix p > 2 and assume that the desired result holds at order p — 1. Use the 
notation u < v for the lexicographical order on hi (in contrast with u -< v for 
the genealogical order !). As usual u < v if u < v and u ^ v. The point is to 
observe that, for every fixed u G the random set 

9 n {v e U : v < u} 

is measurable with respect to the a- field a{K v , v < u). This readily follows from 
the construction of 9. As a consequence, the event 



{u p = u}n {#{9) > P } 
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is measurable with respect to a(K v ,v < u). It is also easy to see that the same 
measurability property holds for the event 

{u p = u}n {#(0) < P }. 

Hence {U p = u} is measurable with respect to a(K v ,v < u). 
Finally, if go, gi, ■ ■ ■ , g p are nonnegative functions on {0, 1, . . .}, 

E[go(K Uo )g 1 (K Ul )...g p (K Up )} 



E * 

u <u 1 <---<u p 

E e 

u Q <Ui<---<u p 



l{Ua=u„,...,u P =u p } 9o(K Uo ) . . . g p (K Up ) 

1 { Uo =u , ■ ■ ■ ,u p =u p } 3o (K Uo )...g p -i (K Up _ 1 ) E [g p (K Up )] 



because K Up is independent of a(K v ,v < u p ), and we use the preceding mea- 
surability property. 

Then E[g p (K Up )\ — n{g p ) does not depend on u p , and taking g p — 1 in the 
preceding formula we see that 

E[g (K Uo )9i{K Ul ) . ..g p (K Up )} = S[ 5o (^o)5i(^c/i) ■ • • 9 P -i{K Up-J] V-{g P )- 
An application of the induction assumption completes the proof. □ 

Corollary 1.6 Let (S n , n > 0) be a random walk on Z with initial value S and 
jump distribution v{k) = p(k + 1) for every k > — 1. Set 

T = inf{n > 1 : S n = -1}. 

Then the Lukasiewicz path of a ^-Galton- Watson tree 9 has the same distribution 
as (So, S\, ... , St)- In particular, #(0) and T have the same distribution. 

This is an immediate consequence of the preceding proposition. 



1.3. Convergence to Brownian motion 

Our goal is to show that the height functions (or contour functions) of Galton- 
Watson trees (resp. of Galton- Watson forests) converge in distribution, modulo 
a suitable rescaling, towards Brownian excursions (resp. reflected Brownian mo- 
tions). 

We fix a critical offspring distribution /j, with finite variance a 2 > 0. Note 
that the criticality means that we now have 

oo 

Let #i, 62, ■ ■ ■ be a sequence of independent /x-Galton- Watson trees. With each 
9i we can associate its height function (hg^n),® < n < #(^) — I)). We then 
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define the height process (H n , n > 0) of the forest by concatenating the functions 
/i0D hg 2 , . . .: 

H n = he i (n-(#(0 1 ) + ■•■ + #&-!))) 

if ) + • • • + #(0;_i) < n < #(6»i) + • ■ • + #(6> 4 ). Clearly, the function 
(H n , n > 0) determines the sequence of trees. To be specific, the "fc-th excursion" 
of H from (more precisely, the values of H between its fc-th zero and the next 
one) is the height function of the fc-th tree in the sequence. 

By combining Corollary II . 61 with Proposition II. 21 we arrive at the following 
result (cf Corollary 2.2 in 

Proposition 1.7 We have for every n > 

H n = #{ke{0,l,...,n-l}:S k = inf Sj). (1) 

k<j<n 

where {S ni n > 0) is a random walk with the distribution described in Corollary 

This is the main ingredient for the proof of the following theorem. By defini- 
tion, a reflected Brownian motion (started at the origin) is the absolute value of 
a standard linear Brownian motion started at the origin. The notation [x] refers 
to the integer part of x. 

Theorem 1.8 Let 81,62, ■■■ be a sequence of independent fi-Galton- Watson 
trees, and let (H n , n > 0) be the associated height process. Then 

(-^ff[ pt] ,i>0) ^ (-7t,*>0) 

where 7 is a reflected Brownian motion. The convergence holds in the sense of 
weak convergence on the Skorokhod space D(R+,R+). 

Let us establish the weak convergence of finite-dimensional marginals in the 
theorem. 

Let S = (S n ,n > 0) be as in Proposition ll.7l Note that the jump distribution 
v has mean and finite variance cr 2 , and thus the random walk S is recurrent. 
We also introduce the notation 

M n = sup S k , I n = inf S k ■ 

0<k<n 0<k<n 

Donsker's invariance theorem gives 

(-^b*],t>0) ^ (vB t ,t>0) (2) 

where B is a standard linear Brownian motion started at the origin. 

For every n > 0, introduce the time-reversed random walk S n defined by 

D k ~ °n J (n— k)+ 
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and note that (S^,0 < k < n) has the same distribution as (S n ,0 < k < n). 
From formula Q), we have 

H n = #{ke {0,1,. ..,n-l}:S k = inf S,} = 

k<j<n 

where for any discrete trajectory u> — (oj(0), w(l), . . .), we have set 
$„(w) = € {1, . . . ,n} : uj{k) = sup uj(j)}. 

0<j<k 



We also set 



K n = $„(5) = #{k e {1, . . . ,n} : S k = M fe }. 



Lemma 1.9 Define a sequence of stopping times Tj, j = 0, 1, . . . inductively by 
setting Tq = and /or every j > 1, 

T 3 = inf{n > 7}_i : 5„ = M n }. 

TTien i/ie random variables — St^x, j = 1,2, . . . are independent and iden- 
tically distributed, with distribution 



P[S Tl =k} = v([k,oc[) , k>0. 

Proof. The fact that the random variables St* — Stj_i j j = 1:2,... are indepen- 
dent and identically distributed is a straightforward consequence of the strong 
Markov property. It remains to compute the distribution of St x ■ 

The invariant measure of the recurrent random walk S is the counting mea- 
sure on Z. By a standard result, if Rq = inf{ri > 1 : S n = 0}, we have for every 
i € Z, 

.Ro-i 



E 



!{S„=i} 



1. 



Notice that T\ < Rq and that the random walk takes positive values on ]Xi, Rq[. 
It easily follows that for every i < 



E 



n=0 



{S„=i} 



Therefore, for any function g : Z 

.Ti-l 



i=0 



(3) 



Then, for any function / : Z 



Z-, 



£[/(S T J] = i?[£l {fc<Tl} /(^ +1 ) Ms fc+1 >o} 

fe=0 
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= ^2 E [ l {k<T 1 }f(S k+1 ) l{S fc+1 >0} 
k=0 

oo oo 

= J2 E [ 1 { k < T i}J2 V ^f( Sk +j) 1 {S k +J>Q} 
fe=0 j=0 
— oo oo 

t=0 3=0 

oo oo 

m— j—rn 

which gives the desired formula. In the third equality we used the Markov prop- 
erty at time k and in the fourth one we applied 0. O 

Note that the distribution of St ± has a finite first moment: 

oo oo . , . 1 \ 2 

E[S Tl ] = E fc K[*> «>)) = E ^2 Ki) - y • 

fe=0 j=o 
The next lemma is the key to the first part of the proof. 
Lemma 1.10 We have 

H n (P) 2 



S n - In n ^°° V 2 ' 

(P) 

where the notation — > means convergence in probability. 
Proof. From our definitions, we have 

M n = ( s n ~ St^ ) = E(^ fc - &■»_,)• (4) 

T fc <n fc=l 

Using Lemma ll .91 and the law of large numbers (note that K n — > oo), we get 

M n (a.s.) ct 2 

— > h[b Tl \ = ~tt- 

J\ n n— too Z 

By replacing S with the time-reversed walk S n we see that for every n, the 
pair {M n ,K n ) has the same distribution as (S n — I n ,H n ). Hence the previous 
convergence entails 

S n ~ In (P) 

H n n—KX 2 ' 

and the lemma follows. □ 
From we have for every choice of < t\ < £2 < • • • < t m , 
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W ~[B tl - inf B s 

~ 0<S<*1 



, a - inf B, 



0<s<t 7 , 



However, a famous theorem of Levy states that the process 



It = B t 



inf B s 

0<s<t 



is a reflected Brownian motion. This completes the proof of the convergence of 
finite-dimensional marginals in Theorem ll.8l 

We will now discuss the functional convergence in Theorem 1 1.81 To this end, 
we will need more precise estimates. We will give details of the argument in the 
case when \i has small exponential moments, that is there exists A > such that 

oo 

^e Afc /j(fc) < oo. 

k=0 

Our approach in that case is inspired from |31j . See |28j for a proof in the general 
case. We first state a lemma. 

Lemma 1.11 Let e € (0, |). We can find e' > and an integer N > 1 such 
that, for every n > N and £ G {0, 1, . . . ,n}, 

2 

P[\Mt - a —K f \ > ni +£ ] < exp(-n £ '). 

We postpone the proof of the lemma and complete the proof of Theorem ll.8l 
We apply Lemma fl . 1 II with e — 1/8. Since for every n the pair (M n ,K n ) has 
the same distribution as (S n — I n , H n ), we get that, for every sufficiently large 
n, and I € {0, 1, . . . , n}, 



P\\Si-h 



-Hi\ > n»] < exp(— n £ ). 



Hence 



P 



sup | Se — h 

0<£<n 



-HA > n8 



< n exp(— n e ). 



Let A > 1 be a fixed integer. We deduce from the preceding bound that, for 
every p sufficiently large, 



sup \S [pt] - I M - —H [pt] | > (Ap) « < Ap exp(-(A P y ). (5) 

0<t<A * J 



A simple application of the Borel-Cantelli lemma gives 



sup 

0<t<A 



Si 



[ P t] - i\pt\ 



cr 2 Hfr t ] 



2 Vp 
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It is now clear that the theorem follows from the convergence 

(-j=(S [pt] - / M ),t > 0) ^ (a(B t - o mf B s ),t > 0), 

which is an immediate consequence of J2J). □ 

We still have to prove Lemma fl. Ill We first state a very simple "moderate 
deviations" lemma for sums of independent random variables. 

Lemma 1.12 Let Yi,Y 2 , ... be a sequence of i.i.d. real random variables. We 
assume that there exists a number A > such that i?[exp(A|Yi|)] < oo, and that 
E[Yi] = 0. Then, for every a > 0, we can choose N sufficiently large so that for 
every n > N and £ 6 {1,2,..., n}, 

P[|Yi + • • • + Y e \ > n^ +a ] < exp(-n Q/2 ). 

Proof. The assumption implies that i?[e AYl ] = 1 + cA 2 + o(A 2 ) as A — * 0, where 
c = ivar(Yi). Hence we can find a constant C such that for every sufficiently 
small A > 0, 

E[e XY ^] <e cx \ 
It follows that, for every sufficiently small A > 0, 



P[Yi +... + Y t > n^ +a ] < e- An5+ ° E[e MYl+ - +Y ^} < e^™ 1 



e 



If n is sufficiently large we can take A = n x / 2 and the desired result follows 
(after also replacing Y t with — Yj). □ 

Let us now prove Lemma fl. Ill We choose a £ (0,e/2) and to simplify nota- 
tion we put m n = [n5+°]. Then, for every I € {0, 1, ... , n}, 

2 2 

P[\M t - U —Kt\ > ni +e ] < P[K e > m n ] + P[\M t - a —K t \ > n^ +e ;K e < m n ). 

(6) 

Recalling (J3J), we have first 

2 

P[\M e - ^K t \ > ni +E ;K e < m n ] 
< P 



< P 



k 

sup \J2((S Tj -S Tj _ 1 )-?-)\>ni + * 

0<fe<m„ , * 

sup |£((S T .-S Tj _J-ff)|> m t +£ 

0<fe<m„ ~[ ^ 



■2 



< m n exp(— m\ 



where the last bound holds for n large by Lemma [1.121 Note that we are as- 
suming that fi has small exponential moments, and the same holds for the law 
of S'tj by Lemma [Ql which allows us to apply Lemma H. 121 
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We still need to bound the first term in the right-hand side of JJJ. Plainly, 
P[K e > m n ] < P[K n > m n ] < P[S Tmn < M n ], 

and so 

P[K t > m n ] < P[S Tmn < ] + P[M n > n^+f ]. 

Applying again Lemma 11.121 we get that for n large, 

P[M n > < n sup P[S e > n^+*] < n exp(-n a/4 ). 

l<t<n 

Finally, 

2 2 

P[S Tmn < ni+^} = P[S Tmn - yrn„ < n^+f - %-m n \ 
2 

and since ST mn — \rn n is the sum of m n i.i.d. centered random variables having 
small exponential moments, we can again apply Lemma 11.121 (or a classical 
large deviations estimate) to get the needed bound. This completes the proof of 
Lemma 11.111 □ 



1. 4- Some applications 

Let us first recall some important properties of linear Brownian motion. Let (3 be 
a standard linear Brownian motion started at 0. Then there exists a continuous 
increasing process L\ = called the local time of f) at such that if N e (t) 

denotes the number of positive excursions of [3 away from with height greater 
than e and completed before time t, one has 

Hin2eN E (t) = L° t 

for every t > 0, a.s. The topological support of the measure dL\ coincides a.s. 
with the zero set {t > : /3 t = 0}. Moreover, the above-mentioned Levy theorem 
can be strengthened in the form 

(B t -B t ,-B t ;t>0) { = ] (|A|,i?CS);*>0) 

where B t = inf < s < t B s . See e.g. jSHI Chapter VI, Theorem VI.2.3. 
Keeping the notation of subsection ll.31 we set for every n > 0, 

A„ = k iff + • • • + #(fl fc _i) < n < #(00 + • • • + #(6 k ) 

in such a way that k is the index of the tree to which the n th -visited vertex 
belongs. 
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The convergence stated in Theorem 1 1 . 81 can now be strengthened in the fol- 
lowing form: 

(^^A^jt^O) ^ (-\(3 t \,aL° t (f3);t>0). (7) 
y/P yP p^oo a 

This is a simple consequence of our arguments: It is easily seen that 

A„ = 1 - inf Sj = 1 - I n . 

On the other hand, we saw that, for every A > 0, 
S\pt] - I[ P t] a 2 -ff[ P t] 



sup 

0<t<^4 



— ► , a.s. 

p~ ->oo 



Vp 2 Vp 

Combining with Donsker's theorem, we get 

(^M^Am^ 0) ^ (^(B t -B t ),-aB t ;t> 0) 

and an application of Levy's theorem in the form recalled above yields J7J). 

We will now apply Q to study the asymptotics of a single Galton- Watson 
tree conditioned to be large. We write h(9) = sup{|w| : v £ 9} for the height 
of the tree 9. Let us fix x > and for every integer p > 1 denote by 0{ x Vp) a 
random tree with distribution 

11^ (da \ h(a) > xVp) 

where we recall that ±1^ is the law of the Galton- Watson tree with offspring 
distribution fi. 

We denote by ijf^VPl the height function of 9^ x ^ . By convention, H„ x ^ — 
if n> #(6>{ a; VP}). 

Corollary 1.13 We have 

(^H};f\t>0) n ( 2 -e°*'\t>0) 
y/p lP * p— >oo a 

where e ax l 2 is a Brownian excursion conditioned to have height greater than 
ax/2. 

The excursion e ax l 2 can be constructed explicitly in the following way. Set 

T = inf{f > : |/3 t | > ax/2} 
G = sup{t <T:/3 t = 0} 
D = m£{t >T:pt = 0}. 

Then we may take 

ax 12 a I 

e* = \P(G+t)AD\- 
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Proof. We rely on l|7)l. From the Skorokhod representation theorem, we may 
for every p > 1 find a process (H^ , ) such that 

[M t ,A t ) t >o - { — H [pt] ,— A [pt ]jt>o 

and 

(i? t (p) ,A^) t > — > (-|/3 t |,aL?(/3)) t > . (8) 

p— too (7 

uniformly on every compact set, a.s. 
Set 

=inf{i > : |ff t (p) | > x} 
G<*> = sup{t < rW : ={)}-p- 1 
= inf{< > : H^ p) = 0}. 

The reason for the term — in the formula for comes from the integer 
part in -^H[ pt y. We want the process H^ p)+t to stay at during the interval 

By construction, (H^j (p)+t ^ ADip) , t > 0) has the distribution of the (rescaled) 
height process of the first tree in the sequence 0\, 62, ■ ■ ■ with height greater than 
Xy/p, which is distributed as 0{ x Vp} . Therefore, 

(H (p) t > () \ (£) (J_ h { x vp} t > ()) 

The corollary will thus follow from (Jj) if we can prove that 

q{p) a - s -, q ; £)(p) a - s -) d 

Using the fact that immediately after time T the process \(3\ hits levels strictly 
larger than crx/2, we easily get from (JSJ that 

p — too 

From this and (JSJ again it follows that 

liminf D {p) > D a.s. 
limsupG (p) >G a.s. 

Let us prove that we have also lim sup £>(p) < D a.s. (the same argument works 
for liminf G^). Let us fix t > 0. From the support property of local time, we 
have 

L\ > L° D a.s. on {D < t}. 

Thanks to ©, we get 

A[ p) > aL° D for p large, a.s. on {D < t}. 
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Now note that L° D = i^. Since converges to T, © also shows that A^, 
converges to aLj,, and it follows that 

A { t p) > A^} p) for p large, a.s. on {D < t}. 

Observing that stays constant on the interval [T^ p \ D^), we conclude that 

t > for p large, a.s. on {D < t}. 

This is enough to prove that limsupD^- 1 < D a.s. □ 

Replacing p by p 2 and taking x = 1, we deduce from the corollary (in fact 
rather from its proof) that 

p^ p — >oo 

where C<t/2 = D — G is the length of excursion e°"/ 2 . Indeed this immediately 
follows from the convergence of - G^ towards D-G and the fact that, 
by construction 

in the notation of the preceding proof. 

Notice that the Laplace transform of the limiting law is known explicitly : 



£[exp(-AC /2 )] = , [ exp(-aV2A/2). 
smh(oV 2A/2) 

This basically follows from the Williams decomposition of Ito's excursion mea- 
sure (Theorem XII. 4. 5 in [35]) and the known formulas for the hitting time of 
a/2 by a three-dimensional Bessel process or a linear Brownian motion started 
at (see e.g. [55]). 

Exercise. Show the convergence in distribution of p^ 1 h(9^) and identify the 
limiting law. 

We will now discuss "occupation measures" . Rather than considering a single 
tree as above, we will be interested in a finite forest whose size will tend to oo 
with p. Precisely, we fix b > 0, and we set 



H p 



H n if A„ < bp, 
if A„ > bp, 



in such a way that H p is the height process for a collection of [bp] independent 
Galton- Watson trees. Then it easily follows from that 

(-K H] ,t>0) ^ (-\(3\ tATb/ „,t>0), (9) 
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where, for every r > 0, 



Indeed, we can write 



r r := inf{f > : L° t > r}. 



-Hi 



-H, 



b 2 *] ~ /V(*Ar< p) )] 



where r fc (p) = ^ inf{n > : A„ > bp} = inf{t > : -A|p2 t ] > b}. Then 
observe from (J7J) that 

((ii/ [p2t] ) t>0 , r b (p) ) m ((-\Pt\)t>o,r b/<7 ) 

p p^oo (J — ' 

and © follows. 

Taking b = 1, we deduce from 10 that, for every x > 0, 



we 



P 



sup > 

l<i<p 



p 



2 1 2 
sup -\f3 t \ > x — 1 - exp( 5—). 



The last equality is a simple consequence of excursion theory for linear Brownian 
motion (see e.g. Chapter XII in jSHI)- Now obviously 



P 



sup h(0i) >px = 1 - (1 - P[h(9) > px\) p 

l<i<p J 



and we recover the classical fact in the theory of branching processes 

P[h{9)>n] ~ * 
n^oo a n 

We now set Zq = p and, for every n > 1 , 
p 

Z P n = G ^ '■ M = "> = #{ fc >°- HP k= "1- 



From Proposition II. 31 we know that (Z^,n > 0) is a Galton- Watson branching 
process with offspring distribution fj,. We can thus apply the classical diffusion 
approximation to this process. 



Theorem 1.14 



(-2L],t>o) ^ (x t ,t>o), 



where the limiting process X is a diffusion process with infinitesimal generator 
i(j 2 a;-^-T, which 

2 dx z ; 

ential eguation 



\o- 2 x-^, which can be obtained as the unique solution of the stochastic differ- 



dX t = o-^Xt dB t 
X = 1. 
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For a proof, see e.g. Theorem 9.1.3 in Ethier and Kurtz ^H]- It is easy to see 
that, for every p > 1, 

• is a martingale, 

• (ZP) 2 — a 2 J2k=o IS a martingale, 

which strongly suggests that the limiting process X is of the form stated in the 
theorem. 

The process X is called Feller's branching diffusion. When a — 2, this is also 
the zero-dimensional squared Bessel process in the terminology of Note 
that X hits in finite time and is absorbed at 0. 

To simplify notation, let us fix /x with a = 2. Let fi,...,f q heq continuous 
functions with compact support from K + into K_|_. As a consequence of (jHJ we 
have 

(p) 

( r h(-H? H] )dt) ^ ( r /2 M\(3 t \)dt) 

^ JO P ' 1< 1 <9 P^°° \ JQ / l<i<q 

On the other hand, 

(p) 

Jo P p t~o P Jo P P lP 1 

By using Theorem II. 141 we see that 

( [ Tl/2 M\Pt\)dt) ( = } ( ^ 'da Ma) X a ) 

In other words, the occupation measure of over the time interval [0, T1/2], 
that is the measure 

/— > [ Tl/2 f(\Pt\)dt 
Jo 

has the same distribution as the measure X a da. We have recovered one of the 
celebrated Ray-Knight theorems for Brownian local times (see e.g. Theorem 
XI.2.3 in ED). 

1.5. Gallon- Watson trees with a fixed progeny 

We can also use Theorem 11.81 to recover a famous result of Aldous concerning 
Galton- Watson trees conditioned to have a large (fixed) number of vertices (see 
|2], Aldous dealt with the contour function rather than the height function, but 
this is more or less equivalent as we will see in the next subsection). We will 
follow an idea of [2H ■ We assume as in the end of subsection 11.31 that fj, has 
a small exponential moment. Our results hold without this assumption, but it 
will simplify the proof. 

For every p > 1 we denote by 6^ a /x-Galton- Watson tree conditioned to 
have #(#) = p. For this to make sense we need P(#(0) = p) > for every 
p > 1, which holds if /z(l) > (in fact, we only need P(#(0) = p) > for 
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p large, which holds under an aperiodicity condition on /i). In the notation of 
subsection 1 1.21 the distribution of 9^ is n M (c?a | #(a) = p). 

We denote by (H/f^oKkKp the height process of 9^ p \ with the convention 

=0. 

We also need to introduce the normalized Brownian excursion (e t )o<t<i- This 
is simply the Brownian excursion conditioned to have length 1 . For instance, we 
may look at the first positive excursion of (3 (away from 0) with length greater 
than 1, write [G, D] for the corresponding time interval, and set 

Et = P{G+t)AD t t >0 

and 

et = 7w=n E(D - G)t ' 

A more intrinsic construction of the normalized Brownian excursion will be 
presented in the next subsection. 

Theorem 1.15 We have 



(^a,0<t<l) n (^e t ,0<i<l) 

tJp lP 1 p^oo a 



This is obviously very similar to our previous results Theorem ll.Hl and Corol- 
lary^^] However, because the present conditioning "degenerates in the limit" 
p — > oo (there is no Brownian excursion with length exactly equal to 1), we 
cannot use the same strategy of proof as in Corollarv ll.131 

Proof. Let (H n , n > 0) be as in Theorem ll.Sl the height process associated with 
a sequence of independent /i-Galton- Watson trees. We may and will assume that 
H is given in terms of the random walk S as in (QJ. 

Denote by T\ the number of vertices of the first tree in the sequence, or 
equivalently 

Ti = inf{n > 1 : H n = 0} = inf{n > : S n = -1}. 

A simple combinatorial argument (consider all circular permutations of the p 
increments of the random walk S over the interval [0,p]) shows that, for every 
P > 1, 

P(T 1 ^p)^^P(S p = -l). 

On the other hand classical results for random walk (see e.g. P9 in Chapter II 
of 00]) give 

lim ^P(S P = -1) ' 



and it follows that 

P(T 1=P ) i=. (10) 

ay / 2np 3 
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Recall from the end of the proof of Theorem 11.81 (see ©) that we can find 
e > so that, for p large enough 



. Hi, 



sup |-M 

o<t<i VP 



2 S[ pt ] 



iM|> p -l/8 



< exp(-p £ ). 



By comparing with (|10() , we see that we have also for p large 
. -ffrrrfi 2 0[ p t] — /i 

~2 



P 



sup 

o<t<i VP 



Wi >p -i/ 



T x =p 



< exp(-p £ ), 



for any e' < e. Since /„ = for < n < Ti, we have also for p large 



P 



Pi 



sup 



\pt] 



2 Si 



o<t<i vp a Vp 



M , > p -i/ 



T x =p 



< exp(-p £ ). 



Now obviously (Pt , < k < p) has the same distribution as [Hk , < k < p) 
under P(- | T± = p). Therefore Theorem ll.lSl is a consequence of the last bound 
and the following lemma which relates the normalized Brownian excursion to 
the random walk excursion with a fixed long duration. 

Lemma 1.16 The distribution of the process (j^=S[pt], < t < 1) under the 
conditional probability P(- | T± — p) converges as p tends to oo to the law of the 
normalized Brownian excursion. 

We omit the proof of this lemma, which can be viewed as a conditional version 
of Donsker's theorem. See Kaigh |24j . 

See also Duquesne ^U] f° r generalisations of Theorem ll.151 

Application. We immediately deduce from Theorem ll . 15l that . for every x > 0, 

ox s 



lim P(h(9) > x^ | #{6)=p) = P( sup e t > — ). 

P^OO 0<t<l i 

There is an explicit (complicated) formula for the right-hand side of 



(11) 



Combinatorial consequences. For several particular choices of /i, the mea- 
sure H fl (da | #(a) = p) coincides with the uniform probability measure on a 
class of "combinatorial trees" with p vertices, and Theorem 11.151 gives informa- 
tion about the proportion of trees in this class that satisfy certain properties. 
To make this more explicit, consider first the case when is the geometric dis- 
tribution with parameter | (fJ-(k) — 2~ fc_1 ). Then U ll (da | #(a) = p) is the 
uniform distribution on the set A p of all rooted ordered trees with p vertices 
(this follows from Proposition II. 4fl . Thus (|llf) shows that the height of a tree 
chosen at random in A p is of order ^fp, and more precisely it gives the asymp- 
totic proportion of those trees in Ap with height greater than x^Jp. Similar 
arguments apply to other functionals than the height: For instance, if < a < b 
are given, we could derive asymptotics for the number of vertices in the tree be- 
tween generations a^fp and b^/p, for the number of vertices at generation [a^/p] 
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that have descendants at generation etc. The limiting distributions are 

obtained in terms of the normalized Brownian excursion, and are not always 
easy to compute explicitly ! 

Similarly, if u is the Poisson distribution, that is n(k) = e _1 /fc!, then Tl^da \ 
#(a) = p) yields the uniform distribution on the set of all rooted Cayley trees 
with p vertices. Let us explain this more in detail. Recall that a Cayley tree 
(with p vertices) is an unordered tree on vertices labelled 1,2, ... ,p. The root 
can then be any of these vertices. By a famous formula due to Cayley, there are 
p p ~ 2 Cayley trees on p vertices, and so p p ~ 1 rooted Cayley trees on p vertices. If 
we start from a rooted ordered tree distributed according to U, 1 (da | #(a) = p), 
then assign labels 1,2, ... ,p uniformly at random to vertices, and finally "forget" 
the ordering of the tree we started from, we get a random tree that is uniformly 
distributed over the set of all rooted Cayley trees with p vertices. Hence Theorem 
11.151 and in particular (|llfl also give information about the properties of large 
Cayley trees. 

As a last example, we can take fi — ^(^0 + ^2)7 and it follows from Proposition 
11.41 that, provided p is odd, 11^ (da | #(a) = p) is the uniform distribution over 
the set of (complete) binary trees with p vertices. Strictly speaking, Theorem 
11.151 does not include this case, since we assumed that /i(l) > 0. It is however 
not hard to check that the convergence of Theorem 1 1 . 1 51 still holds, provided we 
restrict our attention to odd values of p. 

It is maybe unexpected that these different classes of combinatorial trees give 
rise to the same scaling limit, and in particular that the limiting law appearing in 
l|lll) is the same in each case. Note however that the constant a varies: a 2 = 1 for 
(complete) binary trees or for Cayley trees, whereas a 2 = 2 for rooted ordered 
trees. 

As a final remark, let us observe that the convergence in distribution of 
Theorem 1 1.1 51 is often not strong enough to deduce rigorously the desired com- 
binatorial asymptotics (this is the case for instance if one looks at the height 
profile of the tree, that is the number of vertices at every level in the tree). 
Still Theorem 11.151 allows one to guess what the limit should be in terms of 
the normalized Brownian excursion. See in particular [H] for asymptotics of the 
profile that confirmed a conjecture of Aldous. 

1 . 6. Convergence of contour functions 

In this subsection, we briefly explain how the preceding results can be stated 
as well in terms of the contour processes of the trees rather than the height 
processes as discussed above. The contour function of a tree was discussed in 
subsection 11.11 (see Fig0 . Notice that in contrast to the height process it is 
convenient to have the contour function indexed by a real parameter. 

We will give the result corresponding to Theorem 1 1.81 So we consider again a 
sequence 61,62, ■■ ■ of independent /x-Galton- Watson trees and we denote by 
(Ct,t > 0) the process obtained by concatenating the contour functions of 
6\ , 62 , • ■ ■ Here we need to define precisely what we mean by concatenation. 
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In FigQ]and the discussion of subsection ll.il the contour function of a tree is 
naturally defined on the time interval [0, ((6)}, where ((0) = 2(#(0) - 1). This 
has the unpleasant consequence that the contour function of the tree consisting 
only of the root is trivial. For this reason we make the slightly artificial conven- 
tion that the contour function C t {9) is defined for < t < £(0) = 2#(0) — 1, by 
taking C t — if £(0) < t < £(0). We then obtain (Ct,t > 0) by concatenating 
the functions (C t (0i),O < t < £(0i)), (C t (0 2 ),O < t < £(0 2 )), etc. 
For every n > 0, we set 

J„ = In - H„ + I n . 

Note that the sequence J n is strictly increasing and J n > n. 

Recall that the value at time n of the height process corresponds to the 
generation of the individual visited at step n, assuming that individuals are 
visited in lexicographical order one tree after another. It is easily checked by 
induction on n that [J n , J n +i] is exactly the time interval during which the 
contour process goes from the individual n to the individual n + 1. From this 
observation, we get 

sup \C t -H n \ < \H n+1 -H n \ + l. 

te[J n ,Jn+l] 

A more precise argument for this bound follows from the explicit formula for Ct 
in terms of the height process: For t S [J n , J n +i], 

C t = H n - (t - J n ) if t€[J n , J n+ 1 - 1] , 

c t = {H n+1 - (j„ +1 - if * e [J n+1 - 1, J n+1 }. 

These formulas are easily checked by induction on n. 

Define a random function cp : M. + — > {0,1,...} by setting ip(t) = n iff 
t € [J n , Jn+i)- From the previous bound, we get for every integer m > 1, 

sup \C t - H v(t) \ < sup \C t - H v{t) \ < 1 + sup \H n+1 - H n \. (12) 

te[0,m] ie[0,J m ] n<m 

Similarly, it follows from the definition of J n that 

sup \ip{t)-U< sup \(p(t) -U<- sup H n + hl m \ + l. (13) 

te[0,m] ^ t£[0,J m ] ^ ^ n<m ^ 

Theorem 1.17 We have 

-^C 2pt ,t>o) ^> (~|/? t |,*>0). (14) 
where (3 is a standard linear Brownian motion. 

Proof. For every p > 1, set </? p (t) = p^ 1 tp(pt). By (|12|) . we have for every m > 1, 



sup 



1 1 



< -p + -p sup |if[pt]+i - H[pt]| — ► (15) 
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in probability, by Theorem II. 81 

On the other hand, the convergence J3J) implies that, for every m > 1, 



1 j _W 
, — 'mi) 
^/p r p-+oo t< 



Then, we get from 



I mp a mi B t . (16) 

1 I 2 

sup \(p p (2t) -t\<- sup i/ fc + -|/ 2mp | + > (17) 

t<m P fc<2mp P P 

in probability, by Theorem 11.81 and (|l(j|| . 

The statement of the theorem now follows from Theorem 1 1.81 <|15|) and (|17|) . 

□ 

Remark. There is one special case where Theorem 11.171 is easy, without any 
reference to Theorem ll.8l This is the case where /j, is the geometric distribution 
/i(fc) = 2~ k ~ 1 , which satisfies our assumptions with a 2 = 2. In that case, it is not 
hard to see that away from the origin the contour process (C n ,n > 0) behaves 
like simple random walk (indeed, by the properties of the geometric distribution, 
the probability for an individual to have at least n + 1 children knowing that he 
has at least n is 1/2 independently of n) . A simple argument then shows that 
the statement of Theorem 1 1 . 1 71 follows from Donsker's invariance theorem. 

Clearly, Corollarv ll . 1 31 and Theorem l 1 . 1 51 can also be restated in terms of the 
contour process of the respective trees. Simply replace H^) by Wwi (with an 

obvious notation) in Corollary 11.131 and H^l by Cnpt] m Theorem II .151 



1 . 7. Conclusion 



The various results of this section show that the rescaled height processes (or 
contour processes) of large Galton- Watson trees converge in distribution towards 
Brownian excursions. Still we did not assert that the trees themselves converge. 
In fact, a precise mathematical formulation of this fact requires a formal defi- 
nition of what the limiting random trees are and what the convergence means. 
In the next section, we will give a precise definition of continuous trees and dis- 
cuss a topology on the space of continuous trees. This will make it possible to 
reinterpret the results of this section as convergence theorems for random trees. 



Bibiographical notes. The coding of discrete trees by contour functions (Dyck 
paths) or Lukasiewicz words is well known: See e.g. Theorem 1 1 . 81 can be 
viewed as a variant of Aldous' theorem about the scaling limit of the contour 
function of Galton- Watson trees [2] • The method that is presented here is taken 
from |2Hj, with an additional idea from [2]. More general statements can be 
found in Chapter 2 of the monograph See Chapters 5 and 6 of Pitman 
|36|. and the references therein, for more results about the connections between 
trees, random walks and Brownian motion. 
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2. Real Trees and their Coding by Brownian Excursions 

In this section, we first describe the formalism of real trees, which can be used 
to give a precise mathematical meaning to the convergence of rescaled discrete 
trees towards continuous objects. We then show how a real tree can be coded 
by a continuous function in a way similar to the coding of discrete trees by their 
contour functions. Aldous' Continuum Random Tree (the CRT) can be defined 
as the random real tree coded by a normalized Brownian excursion. For every 
integer p > 1, we then compute the p-dimensional marginal distribution (that 
is the law of the reduced tree consisting of the ancestral lines of p individuals 
chosen uniformly at random) of the tree coded by a Brownian excursion under 
the Ito excursion measure. Via a conditioning argument, this leads to a simple 
derivation of the marginal distributions of the CRT. 

2.1. Real trees 

We start with a formal definition. In the present work, we consider only compact 
real trees, and so we include this compactness property in the definition. 

Definition 2.1 A compact metric space (T, d) is a real tree if the following two 
properties hold for every a, 6 G T . 

(i) There is a unique isometric map f a ,b from [0, d(a, b)] into T such that 
/a,fc(0) = a and f a ,b(d(a,b)) = b. 

(ii) If q is a continuous infective map from [0,1] into T, such that q(0) = a 
and q(l) = b, we have 

g([0,l]) = / Oj6 ([0 1 d(a > 6)]). 

A rooted real tree is a real tree (T, d) with a distinguished vertex p = p(T) 
called the root. In what follows, real trees will always be rooted, even if this is 
not mentioned explicitly. 

Let us consider a rooted real tree (T,d). The range of the mapping f a j, in 
(i) is denoted by [[a, b]] (this is the line segment between a and b in the tree) . In 
particular, [[/?, a]] is the path going from the root to a, which we will interpret 
as the ancestral line of vertex a. More precisely we define a partial order on the 
tree by setting a =4 b (a is an ancestor of b) if and only if a G [[p, b]]. 

If a, b G T, there is a unique c G T such that [[p, a]] n [[p, b]] — [[p, c]]. We write 
c = a A b and call c the most recent common ancestor to a and b. 

By definition, the multiplicity of a vertex a G T is the number of connected 
components of T\{a}. Vertices of T\{p} which have multiplicity 1 are called 
leaves. 

Our goal is to study the convergence of random real trees. To this end, it is 
of course necessary to have a notion of distance between two real trees. We will 
use the Gromov-Hausdorff distance between compact metric spaces, which has 
been introduced by Gromov (see e.g. ^1) m view of geometric applications. 
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If (E,S) is a metric space, we use the notation SHaus(K, K') for the usual 
Hausdorff metric between compact subsets of E : 



'Haus 



(K, K') = inf{e > : K C U e (K') and K' C U £ (K)}, 



where U £ (K) :={ie£: 6(x, K) < e}. 

Then, if T and T' are two rooted compact metric spaces, with respective 
roots p and p' , we define the distance cIgh(1~,T') by 

d Gff (T, T') = m£{S Haus (cp(T), ip'(T')) V %(/>), v'(p'))} 

where the infimum is over all choices of a metric space (E, S) and all isometric 
embeddings <p : T — ► E and tpf : T' — > E of T and T into J). 

Two rooted compact metric spaces T\ and are called equivalent if there 
is a root-preserving isometry that maps 1\ onto T^- Obviously <1gh{T~, T') only 
depends on the equivalence classes of T and T'. Then cIqh defines a metric on 
the set of all equivalent classes of rooted compact metric spaces (cf arj d 
(17)1. We denote by T the set of all (equivalence classes of) rooted real trees. 

Theorem 2.1 The metric space (T,cIgh) is complete and separable. 

We will not really use this theorem (see however the remarks after Lemma 
12 .4(1 . So we refer the reader to J7|, Theorem 1 for a detailed proof. 

We will use the following alternative definition of oIgh- First recall that if 
(Ti,d%) and (T2, da) are two compact metric spaces, a correspondence between 
Ti and 72 is a subset TZ of T\ x T2 such that for every xi G T\ there exists 
at least one xi G such that (^1,^2) G 7?. and conversely for every 1/2 € ^2 
there exists at least one y\ G T\ such that (2/1,2/2) G 7?.. The distortion of the 
correspondence TZ is defined by 

dis(7£) = sup{|di(xi,yi) - d 2 (x 2 ,V2)\ ■ (xi,x 2 ), (2/1,2/2) G TZ}. 

Then, if T and T' are two rooted compact metric spaces with respective roots 
p and p' , we have 

d GH (T,T') = l inf dis(ft), (18) 

where C(T, T') denotes the set of all correspondences between T and T (sec 
Lemma 2.3 in |17| - actually this lemma is stated for trees but the proof applies 
as well to compact metric spaces). 



2.2. Coding real trees 

In this subsection, we describe a method for constructing real trees, which is 
particularly well-suited to our forthcoming applications to random trees. We 
consider a (deterministic) continuous function g : [0, 00 [ — ► [0, 00 [ with compact 
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support and such that g(0) = 0. To avoid trivialities, we will also assume that 
g is not identically zero. For every s, t > 0, we set 

m(s,t) = inf g(r), 

re[sAt,sVt] 

and 

d g( s > *) = 5( s ) + 9(t) - 2m g (s, t). 
Clearly d g (s, t) = d g (t, s) and it is also easy to verify the triangle inequality 

dg(s,U) < dg(s,t) +dg(t,U) 

for every s, t, u > 0. We then introduce the equivalence relation s ~ t iff 
d g (s,t) = (or equivalently iff g(s) = g(t) = m g (s,t)). Let T g be the quo- 
tient space 

r 5 = [o,oo[/~. 

Obviously the function d g induces a distance on T g , and we keep the notation 
d g for this distance. We denote by p g : [0, oo[ — > T g the canonical projection. 
Clearly p g is continuous (when [0, oo[ is equipped with the Euclidean metric and 
T g with the metric d g ). 

We set p = p g (Q). If C > is the supremum of the support of g, we have 
p g (t) = p for every t > Q. In particular, T g = p g ([0, Q) is compact. 

Theorem 2.2 The metric space (T g ,d g ) is a real tree. We will view (T g ,d g ) as 
a rooted tree with root p = p g (0)- 

Remark. It is also possible to prove that any (rooted) real tree can be repre- 
sented in the form T g . We will leave this as an exercise for the reader. 



9(r) 




Figure 2 



To get an intuitive understanding of Theorem 12.21 the reader should have 
a look at Figure El This figure shows how to construct a simple subtree of 
T g , namely the "reduced tree" consisting of the union of the ancestral lines 
in Tg of three vertices p g (s),p g (t),p g (u) corresponding to three (given) times 
s,t, u € [0, (}. This reduced tree is the union of the five bold line segments that 
are constructed from the graph of g in the way explained on the left part of 
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the figure. Notice that the lengths of the horizontal dotted lines play no role in 
the construction, and that the reduced tree should be viewed as pictured on the 
right part of Figure El The ancestral line of p g {s) (resp. p g (t),p g (u)) is a line 
segment of length g(s) (resp. g(t),g(u)). The ancestral lines of p g (s) and p g (t) 
share a common part, which has length m g (s, t) (the line segment at the bottom 
in the left or the right part of Figure EJ, and of course a similar property holds 
for the ancestral lines of p g (s) and p g {u), or of p g (t) and p g {u). 

We present below an elementary proof of Theorem 12.21 which uses only the 
definition of a real tree, and also helps to understand the notions of ancestral 
line and most recent common ancestor in T g . Another argument depending 
on Theorem 12.11 is presented at the end of the subsection. See also ^Hj for a 
short proof using the characterization of real trees via the so-called four-point 
condition. 

Before proceeding to the proof of the theorem, we state and prove the fol- 
lowing root change lemma, which is of independent interest. 

Lemma 2.3 Let sq E [0, ([. For any real r > 0, denote by r the unique element 
of [0, £[ such that r — r is an integer multiple of £. Set 

g'(s) = g(s a ) + g(s + s) - 2m g (s , s + s), 

for every s G [0,£], and g'(s) — for s > £. Then, the function g' is contin- 
uous with compact support and satisfies g'(0) — 0, so that we can define T g > . 
Furthermore, for every s,tE [0, £], we have 

d g >(s,t) =dg(so + s,s + t) (19) 

and there exists a unique isometry R from T g i onto T g such that, for every 

«e[o,C], ' 

R(Pg'(s))=Pg(so + s). (20) 

Assuming that Theorem l2.2l is proved, we see that T g > coincides with the real 
tree T g re- rooted at p g (so). Thus the lemma tells us which function codes the 
tree T g re-rooted at an arbitrary vertex. 

Proof. It is immediately checked that g' satisfies the same assumptions as g, so 
that we can make sense of T g i . Then the key step is to verify the relation (|19|) . 
Consider first the case where s, t G [0, C — So[. Then two possibilities may occur. 

If m g (s a + s, s Q + t)> m g (s , s + s), then m g (s a , s + r) = m g (s , s + s) = 
w s (so, so + t) for every r € [s, i], and so 

m g '(s,t) = g(s Q ) + m g (s + s,s +t)- 2m g (s , s + s). 

It follows that 

d g i(s,t) = g'(s) + g'(t) - 2m g >(s,t) 

= g(so + s) - 2m g (s Q , s + s) + g(s + 1) - 2m g (s ,s + t) 

-2(m g (s Q + s, s + 1) - 2m g (s , s + s)) 
= g(s + s) + g(s + 1) - 2m g (s + s, s + t) 
= dg(so+s,s +t). 
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If m g (so + s, so + i) < m g (so, So + s), then the minimum in the definition of 
m g '(s,t) is attained at r\ defined as the first r G [s,t] such that g(so + r) = 
w g (so, so + s) (because for r G [r%, t] we will have g(so + r) — 2m 9 (so, so + r) > 
—m g (sa, s + r) > —m g (so, sq + r*i)). Therefore, 

m g >(s,t) = g(s Q ) - m g (s Q , s Q + s), 

and 

d g t(s,t) = g(so + s)-2m g (s ,s + s)+g(so+t) 
-2ra g (s , s + t) + 2m g (s , s + s) 
= d 5 (s + s,s + t). 

The other cases are treated in a similar way and are left to the reader. 

By lpn?)l. if s, t € [0, £] are such that d g >(s, t) = 0, we have d g (so + s, sq + t) — 
so that p g (so + s) — p g {so + t). Noting that T g i = p g '([0, C]) (the supremum 
of the support of g' is less than or equal to £) , we can define R in a unique way 
by the relation J5DJ). From (|19|l . R is an isometry, and it is also immediate that 
R takes T g * onto T g . □ 

Proof of Theorem 12.21 Let us start with some preliminaries. For a, a' G T g , 

we set er =^ a 1 if and only if d g (a,a') = d g {p,a') — d g {p,a). If er = p 9 (s) and 
cr' = p g (t), it follows from our definitions that er =<! cr' iff m g (s,t) = g(s). It is 
immediate to verify that this defines a partial order on T g . 
For any ao , a G 7^ , we set 

[[<T , cr]] = Ke^ 9 : d fl ((T , cr) = <i g ((To, c') + c? s ((t', a)}. 

If cr = p s (s) and cr' = p g (t), then it is easy to verify that [[p, cr]] fl [[p, cr']] = [[p, 7]], 
where 7 = p g (r), if r is any time which achieves the minimum of g between s 
and t. We then put 7 = cr A cr'. 

We set T 5 [cr] := {cr' e T g : cr ^ cr'}. If T g [cr] 7^ {cr} and cr 7^ p, then T g \T g [a] 
and 7g[cr]\{cr} are two nonempty disjoint open sets. To see that T g \T g [a] is open, 
let s be such that p g (s) = cr and note that T g [a] is the image under p g of the 
compact set {u e [0, £] : m ff (s, u) = (?(s)}. The set 7^[cr]\{cr} is open because if 
cr' G ^[cr] and cr' 7^ cr, it easily follows from our definitions that the open ball 
centered at a' with radius d g (cr, cr') is contained in 7 s [cr]\{cr}. 

We now prove property (i) of the definition of a real tree. So we fix cri and 02 
in T g and we have to prove existence and uniqueness of the mapping f ai ^ . By 
using Lemma \2 . 31 with s such that p g (s ) = o\, we may assume that a\ = p. If 
cr G T g is fixed, we have to prove that there exists a unique isometric mapping 
/ = fp.a from [0, d g (p, cr)] into 7^ such that /(0) = p and f{d g {p,a)) = a. Let 
s G p^ 1 ({cr}), so that g(s) — d g (p, a). Then, for every a G [0, d g (p, cr)], we set 

v(a) = inf{r G [0, s] : m g (r, s) = a}. 

Note that g(v(aj) — a. We put f(a) = p g (v(a)). We have /(0) = p and 
f{d g {p,a)) — cr, the latter because m g (v{g{s)), s) — g(s) implies P g (v(g(s))) = 
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p g (s) = a. It is also easy to verify that / is an isometry: If a, b G [0, d g (p, cr)] 
with a < 6, it is immediate that m g (v(a), v(b)) — a, and so 

d fl (/(o), /(&)) = ff(«(o)) + g(v(b)) -2a = b-a. 

To get uniqueness, suppose that / is an isometric mapping satisfying the 
same properties as /. Then, if a G [0, d g (p, cr)], 

dg{f{ a )> a ) = d g(P> o)-a = d g (p, a) - d g (p, f{a)). 

Therefore, /(a) cr. Recall that a = p g (s), and choose t such that p g (t) = /(a). 
Note that g(t) = d g {p,p g {i)) — a. Since /(a) =<! a we have g(t) = m g (t,s). On 
the other hand, we also know that a — g(v(a)) — m g (y(a), s). It follows that 
we have a — g(t) — g(v(a)) — m g (v(a),t) and thus d g (t,v(a)) = 0, so that 
/(a) = p g (t) = p g (v(aj) = /(a). This completes the proof of (i). 

As a by-product of the preceding argument, we see that /([0, d g (p, a)]) = 
[[p, ct]]: Indeed, we have seen that for every a G [0,d g (p,a)], we have f(a) =4 a 
and, on the other hand, if 7/ =<! a, the end of the proof of (i) just shows that 
V = f(d g (p, ??))■ 

We turn to the proof of (ii). We let q be a continuous injective mapping from 
[0, 1] into T g , and we aim at proving that q([0, 1]) = / g ( ), g (i)([0, d g (q(Q),q(l))]). 
From Lemma |2.rtl again, we may assume that g(0) = p, and we set a = q(l). 
Then we have just noticed that /o, CT ([0, d g (p, a)]) — [[p,cr]]. 

We first argue by contradiction to prove that [[p, a]] C g([0, 1]). Suppose that 
r) € [[p, cr]]\g([0, 1]), and in particular, r\ ^ p.a. Then g([0, 1]) is contained in 
the union of the two disjoint open sets T g \T g [rj\ and T g [ri]\{vi\, with g(0) = p G 
T g [rf\\{ri} and g(l) = a E T g \T g [i]]. This contradicts the fact that g([0,l]) is 
connected. 

Conversely, suppose that there exists a £ (0,1) such that q(a) £ [[p, cr]]. Set 
rj = q(a) and let 7 = a A rj. Note that 7 € [[/?, 77]] n [[77, cr]] (from the definition 
of a A T], it is immediate to verify that d g {rj,a) = d g (rj,j) + ^(7, cr)). From 
the first part of the proof of (ii), 7 € g([0, a]) and, via a root change argument, 
7 G g([a, 1]). Since o is injective, this is only possible if 7 — q(a) = 77, which 
contradicts the fact that 77 ^ [[p, cr]] . □ 

Once we know that (T g , d g ) is a real tree, it is straightforward to verify that 
the notation a =^ cr', [[cr, cr']], crAcr' introduced in the preceding proof is consistent 
with the definitions of subsection 12.11 stated for a general real tree. 

Let us briefly discuss multiplicities of vertices in the tree T g . If a G T g is not 
a leaf then we must have i{a) < r(cr), where 

%) := inf ^({cr}) , r(cr) := sup ^({cr}) 

are respectively the smallest and the largest element in the equivalence class 
of cr in [0,C]- Note that m g (£(a),r((r)) = g(£(a)) = g{r{a)) — d g (p,a). Denote 
by (cii,bi),i G 1 the connected components of the open set (£(a),r(a)) n {t G 
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[0, oo[: g(t) > d g (p, cr)} (the index set X is empty if a is a leaf). Then we claim 
that the connected components of the open set 7^\{<t} are the sets p g ((a,i, h)), 
i 6 1 and T g \T g [a] (the latter only if a is not the root). We have already 
noticed that T g \T g [a] is open, and the argument used above for 7g[cr]\{cr} also 
shows that the sets p g ((a,i,bi)), i £ 1 are open. Finally the sets p g ((a.i,bi)) 
are connected as continuous images of intervals, and T g \T g [a] is also connected 
because if a', a" S T g \T g [cr], [fp, cr']} U [fp, a"]} is a connected closed set contained 
in T g \T g [a}. 

We conclude this subsection with a lemma comparing the trees coded by two 
different functions g and g' . 

Lemma 2.4 Let g and g' be two continuous functions with compact support 
from [0, oof into [0, oof, such that g(0) — g'(0) = 0. Then, 

d GH (T g ,T g >)<2\\g-g% 

where \\g — g'\\ stands for the uniform norm of g — g' . 

Proof. We rely on formula (|18|) for the Gromov-Hausdorff distance. We can 
construct a correspondence between T g and T g i by setting 

TZ = {(u, cr') : cr = p g (t) and a' — p g >(t) for some t > 0}. 

In order to bound the distortion of 7Z, let (cr, a') £ TZ and (rj, r/') G 1Z. By 
our definition of TZ we can find s,t > such that p g (s) = a, p g >{s) = a' and 
p g {t) — rj, p g '{t) — rj . Now recall that 

dgfav) = 9(s)+g(t)-2m g (s,t), 
d gl (a',ri') = g'(s)+g'(t)-2m g ,( Sl t), 

so that 

\d g (a,r))-d g ,(o-', V ')\<4\\g-g'\\. 

Thus we have dis(7?.) < 4\\g — g'\\ and the desired result follows from 1)18(1 . □ 

Lemma 12.41 suggests the following alternative proof of Theorem 12.21 Denote 
by Coo the set of all functions g : [0, oof — ► [0, oof that satisfy the assumptions 
stated at the beginning of this subsection, and such that the following holds: 
There exist e > and p > such that, for every i £ N, the function g is linear 
with slope p or — p over the interval [(i — l)e, is)]. Then it is easy to see that T g 
is a real tree if g £ Coo- Indeed, up to a simple time-space rescaling, g will be 
the contour function of a discrete tree t G A, and T g coincides (up to rescaling) 
with the real tree that can be constructed from t in an obvious way. Then, a 
general function g can be written as the uniform limit of a sequence (g n ) in Coo, 
and Lemma 12.41 implies that T g is the limit in the Gromov-Hausdorff distance 
of the sequence T 9n . Since each T 9n is a real tree, T g also must be a real tree, by 
Theorem 12 . 1 1 (we do not really need Theorem 12 .11 but only the fact that the set 
of all real trees is closed in the set of all compact metric spaces equipped with 
the Gromov-Hausdorff metric, cf Lemma 2.1 in |17j). 
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Remark. Recalling that any rooted real tree can be represented in the form T g , 
the separability in Theorem l2.1l can be obtained as a consequence of Lemma EOl 
and the separability of the space of continuous functions with compact support 
on K+. 

2.3. The continuum random tree 

We recall from Section ^ the notation e = (e t ,0 < t < 1) for the normalized 
Brownian excursion. By convention, we take e t = if t > 1. 

Definition 2.2 The continuum random tree (CRT) is the random real tree T e 
coded by the normalized Brownian excursion. 

The CRT T e is thus a random variable taking values in the set T. Note that 
the measurability of this random variable follows from Lemma 12.41 

Remark. Aldous PQ,|2] uses a different method to define the CRT. The preced- 
ing definition then corresponds to Corollary 22 in [5] . Note that our normaliza- 
tion differs by an unimportant scaling factor 2 from the one in Aldous' papers: 
The CRT there is the tree Te instead of T e . 

We can restate many of the results of Section in terms of weak convergence 
in the space T. Rather than doing this in an exhaustive manner, we will give a 
typical example showing that the CRT is the limit of rescaled "combinatorial" 
trees. 

Recall from subsection 11.11 the notation A for the set of all (finite) rooted 
ordered trees, and denote by A„ the subset of A consisting of trees with n 
vertices. We may and will view each element t of A as a rooted real tree: 
Simply view t as a union of line segments of length 1 in the plane, in the way 
represented in the left part of Figure^ equipped with the obvious distance (the 
distance between a and a' is the length of the shortest path from a to a' in the 
tree). Alternatively, if (Ct,t > 0) is the contour function of the tree, this means 
that we identify t = Tq (this is not really an identification, because the tree t 
has an order structure which disappears when we consider it as a real tree). 

For any A > and a tree T 6 T, the tree AT is the "same" tree with all 
distances multiplied by the factor A (if the tree is embedded in the plane as 
suggested above, this corresponds to replacing the set T by AT). 

Theorem 2.5 For every n > 1, let T(„) be a random tree distributed uniformly 
over A„. Then (2n) _1 / 2 T(„) converges in distribution to the CRT T e , in the 
space T. 

Proof. Let 9 be a Galton- Watson tree with geometric offspring distribution 
fi(k) = 2~ fe_1 , and for every n > 1 let 8 n be distributed as 8 conditioned to 
have n vertices. Then it is easy to verify that 6 n has the same distribution as 
T( n ). On the other hand, let (C™,f > 0) be the contour function of 9 n , and let 



Cr = (2n)- 1 /2 C ' 2 " nt , t>0. 
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From Theorem 11.151 (restated in terms of the contour function as explained in 
subsection I l.(jf> . we have 

(C t V>0) ^ (e t ,i>0). 

n — >oo 

On the other hand, from the observations preceding the theorem, and the fact 
that 9 n has the same distribution as T( n ) i it is immediate that the tree 
coded by C n has the same distribution as (2n) _1 / 2 T(„). The statement of the 
theorem now follows from the previous convergence and Lemma 12.41 □ 

We could state analogues of Theorem 12.51 for several other classes of com- 
binatorial trees, such as the ones considered at the end of subsection 11.51 For 
instance, if t„ is distributed uniformly among all rooted Cayley trees with n 
vertices , then (4n)~ 1 / 2 r„ converges in distribution to the CRT T e , in the space 
T. Notice that Cayley trees are not ordered, but that they can be obtained from 
(ordered) Galton- Watson trees with Poisson offspring distribution by "forget- 
ting" the order, as was explained at the end of subsection ll.5l By applying the 
same argument as in the preceding proof to these (conditioned) Galton- Watson 
trees, we get the desired convergence for rescaled Cayley trees. 

2. 4- The ltd excursion measure 

Our goal is to derive certain explicit distributions for the CRT, and more specif- 
ically its so-called finite-dimensional marginal distributions. For these calcula- 
tions, we will need some basic properties of Brownian excursions. Before dealing 
with the normalized Brownian excursion, we will consider Ito's measure. 

We denote by (B t ,t > 0) a linear Brownian motion, which starts at x under 
the probability measure P x . We set 

S t = sup B s , I t = inf B s 

s<t s <* 

and, for every a € K, T a = inf{t > : B t = a}. The reflection principle gives 
the law of the pair (S t , B t ): If a > and b e (— oo, a], 

P [S t >a,B t <b}= P [B t > 2a - b}. 

It follows that, for every t > 0, the density under Pq of the law of the pair 
(S u B t ) is 

. 2(2o-6) / (2a-6) 2 \ 

7 t (a, b) = exp ( - k_^L) 1 {0 > , 6 < 0} . (21) 

The reflection principle also implies that St and \B t \ have the same distribution. 
Let a > 0. Observing that {T a < t} = {St > a}, Pq a.s., we obtain that the 
density of T a under Pq is the function 

a . a 2 

qa{t) = eM ~^- 
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Notice the relation 7* (a, b) = 2 q2 a -b(t) for a > and b < a. 

For every e > 0, denote by v £ the law of the first excursion of B away 
from that hits level s. More specifically, if G £ = sup{t < T £ : B t = 0} and 
D £ = inf{t > T £ : B t = 0}, v £ is the law of {B^ Ge+t ^ DEl t > 0). The measure 
v £ is thus a probability measure on the set C = C(R.+ ,ffi.+) of all continuous 
functions from R + into R+, and is supported on C £ = {e € C : supe(s) > e}. If 
< e < e', we have 

u e (C e >)=P e [T e ,<T ] = 



and 



For every e > 0, set 



= i/ e (- |C e /) = -M'HC e 



1 

2e 



Then n e / = n e (-nC e /) for every < e < e' . This leads to the following definition. 
Definition 2.3 The a-finite measure n on C defined by 

n = lim T n E 

is called the ltd measure of positive excursions of linear Brownian motion. 

Let us briefly state some simple properties of the Ito measure. First n(de) is 
supported on the set £ consisting of all elements e£C which have the property 
that there exists a = a(e) > such that e(t) > if and only if < t < 
a (the number cr(e) is called the length, or the duration of excursion e). By 
construction, n e is the restriction of n to C £ , and in particular n(C £ ) = (2e) _1 . 
Finally, if T £ (e) = inf{* > : e(t) = e}, the law of (e(T e (e) + t),t> 0) under 
n(- I T £ < 00) = v £ is the law of (B t/ ^T ,t > 0) under P £ . The last property 
follows from the construction of the measure v £ and the strong Markov property 
of Brownian motion at time T £ . 

Proposition 2.6 (i) For every t > 0, and every measurable function g : M + — ► 
K+ such that g(0) = 0, 



n{de) g(e(t)) = J dxq x {t)g{x). (22) 

In particular, n(a > t) — n(e(t) > 0) = (27rf)~ 1//2 < 00. Moreover, 

n( dtg(e(t)}) - jf dxg(x). (23) 

(ii) Let t > and let $ cmd ^/ be two nonnegative measurable functions defined 
respectively on C([0,t],R + ) and C. Then, 



1 



n(de) $(e(r), < r < t)*(e(t + r), r > 0) 
= J n(de) $(e(r),0 < r < t) E e{t) [$(B r ^ To ,r > 0)]. 
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Proof, (i) We may assume that g is bounded and continuous and that there 
exists a > such that g{x) — if x < a. Then, by dominated convergence, 



n{de)g(e{t)) = hm / n(de) g(T £ (e) + t) l {Te ( e )<oo} = Urn— E e [g(B tATo )], 

using the property stated just before the proposition. From formula (|21|l we get 

E s \g(B tATo )]=EMBt)l {t <T }] = E [g(e-B t )l {St<e} ] 

f a 

da I dbg(e — b) j t (a,b). 



The first assertion in (i) now follows, observing that q x (t) = 57t(0, — x). The 
identity n(e(t) > 0) = (27rf)~ 1//2 < oo is obtained by taking g(x) = l{:r>o}- The 
last assertion in (i) follows from (|22|l . recalling that the function t — > q x (t) is a 
probability density. 

(ii) We may assume that $ and ^ are bounded and continuous and that there 
exists a <E (Q,t) such that $(w(r),0 < r < t) = if a;(a) = 0. The proof then 
reduces to an application of the Markov property of Brownian motion stopped 
at time T , by writing 



n(de) $(e(r), < r < i)*(e(f + r), r > 0) 

lim | f»(de) l { T e(e) <oo} $(e(T e (e) + r), < r < t)*(e(T s (e) + t + r), r > 0) 
r $(5 rATo) < r < t) *(B (t+s)ATo , s > 0) 
<f (BrAT , < r < t) E BtATo [V(B sATo ,s > 0)] 
!im/ n(de) l { r e ( e )<oc} *(e(T £ (e) + r), < r < t)E e{Te(e)+t) {*(B sATo ,s> 0)] 
(tie) *(e(r), < r < t) E e{t) [f(B rATo ,r > 0)] . 



lim — E t 

e^o 2e 

lim — E, 



In the first and in the last equality, dominated convergence is justified by the 
fact that l{ Te ( e ) <00 }$(e(T £ (e) + r), < r < t) = if <r(e) < a, and the property 
n(a > a) < oo. □ 

2.5. Finite- dimensional marginals under the ltd measure 

If (T, g?) is a real tree with root p, and if xi, . . . , x p G T, the subtree spanned 
by Xi , . . . , x p is simply the set 



T(xi, ...,x p ) = \j[[p,Xi]]. 
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It is easy to see that T(x%, . . . ,x p ), equipped with the distance d, is again a 
real tree, which has a discrete structure: More precisely, T(x\, . . . ,x p ) can be 
represented by a discrete skeleton (which is a discrete rooted tree with p labelled 
leaves) and the collection, indexed by vertices of the skeleton, of lengths of 
"branches" . 

Rather than giving formal definitions for a general real tree, we will concen- 
trate on the case of the tree T g coded by g in the sense of subsection 12.21 

Recall that C denotes the set of all continuous functions from R + into R+ . We 
consider a general continuous function g G C (in contrast with subsection 12.21 
we do not assume that g(0) = and g has compact support). If < t\ < t 2 < 
■ ■ ■ < t p , we will define a "marked tree" 

9(g;ti, ...,t p )= (r(g;ti, . . . ,t p ), (^)«ETO;ti,...,t p )) 

where r(g; t\, . . . , t p ) € A (A is the set of all rooted ordered trees as in Section^) 
and h v > for every v € r(g;tx, . . . ,t p ). An example with p = 3 is given 
in Figure 01 below. In this example, T(g;ti,t2,ts) = {0, 1, 2, (2, 1), (2, 2)} as 
pictured in the right part of Figure|21 and the numbers h v , v S r(g; t\,t2, £3) are 
the lengths of the bold segments as indicated on the left part of the figure. 



g( r ) 




h t 2 t 3 r T(g;ti,t2,t3) 

Figure 3 



To give a precise definition of 9(g; t±, . . . , t p ), we proceed by induction on p. 
lfp = 1, r(g;ti) = {0} and h (g;t 1 ) = g{h). 

Let p > 2 and suppose that the "marked tree" 9(g; t\, . . . ,tj) has been con- 
structed up to order p — 1. Then there exists an integer k G {1, . . . ,p — 1} and 
k integers 1 < i\ < 12 < ■ ■ ■ < ik < P — 1 such that m g (ti, ti+i) = m g (ti,t p ) 
iff i G Set iq = and ik+i — p by convention. For every t G 

{1, 2, . . . , k + 1}, define g G C by the formula 

9*(t) = 3((*V U e _ 1+1 ) At u ) - m g (tx,tp). 

We then let r(g; t%, . . . , t p ) be the concatenation of the trees r(g e ; t^^+i, . . . , ti t ) 
for 1 < I < k + 1: Precisely, 

fc+i 

r(ff;ti,...,t p ) = {0}U |J{£u : u G T(g e ;t ie _ 1+1 , . . . ,t ie )}. 



J.F. Le Gall/Random trees 



282 



Furthermore, if for £ £ {1, . . . , k + 1} 

we define the marks (^)t,e T (g;ti,...,t p ) by setting 

/i„ = h l u , if u = iu , »e T(g e ;t ie _ 1+1 , ...,t ie ), 

and /10 = m g (ti,t p ) 

This completes the construction of the tree by induction. Note that k + 1 is 
the number of children of in the tree 6{g; t\ , . . . , t p ), and m(ti , t p ) is the mark 
of 0. 

If now g satisfies the conditions in subsection 12.21 it is easy to see that 
6(g; ti, ■ ■ ■ , tp) corresponds to the tree T g (p g (ti), . . . ,p g (t p )) spanned by the 
vertices p g (t\), . . . ,p g {t p ) in the tree T g . More precisely, if we attach to ev- 
ery v £ r(g;ti, ... ,t p ) a line segment in the plane with length h v , in such a 
way that the line segments attached to v and to its children share a common 
end (the same for all children of v) and that the line segments otherwise do not 
intersect, the union of the resulting line segments will give a representative of 
the equivalence class of T g (p g (ti), . . . ,p g (t p )). (Note that the order structure of 
t(<?; t\, . . . , t p ) plays no role in this construction.) 

We let A( p ) be the set of all rooted ordered trees with p leaves (a leaf of a 
tree r G A is a vertex u € r with no child, i.e. such that k u (r) = 0, with the 
notation of Section We denote by Q( p ) the set of all marked trees with p 
leaves: Elements of 0( p ) are of the form 9 — (r, (h v ) ver ) where r £ A( p ) and 
h v > for every v e r. The set 6( p ) is equipped with the obvious topology and 
the associated Borel <r-field. We also consider the corresponding sets of binary 
trees: Aj£? is the set of all binary rooted trees with p leaves (and hence 2p — 1 
vertices), and 0^" is the set of marked trees = (r, (h v ) ver ) whose skeleton r 
belongs to A]?'?. Recall that 



«Ag)- (2 "- 2 » ! 



is the Catalan number of order p — 1 . 

Theorem 2.7 The law of the tree 9(e; t\, . . . , t p ) under the measure 

Tl{de) l^KtxK - KtpKcrie^dti . . . dt p 

is 2 p ~ 1 A p , where A p is the uniform measure on defined by 

[ A p (d9)F(9)= fl[dh v F(T,{h v ,veT}). 

It should be clear from our construction that the tree 9(e;t±, . . . ,t p ) only 
depends on the values of e(ii), . . . , e(t p ) and of the successive minima m e (ti, £2), 
m e (t2, £3), ■ ■ ■ , m e (t p _i,t p ). The key tool in the proof of Theorem 12 . 71 will thus 
be the following proposition. To simplify notation we write m(s,t) = m e (s,t). 
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Proposition 2.8 Let F be nonnegative and measurable on R 2 ^ 5-1 . Then 
n( / dti...dtpf(rn(h,t2),...,m(t p - 1 ,t p ),e(t 1 ),...,e(t p ))) 

V J{0<ti< — <t p <ff} ' 

= 2 P ~ 1 / dai..da p -id(3i.. d/3 p ( TT l{ aj <3 i Afl i+1 } ..,a p -x,Px, .. ,f3 p ). 

H*' 1 i="i y 

Before proving this proposition, we state a lemma which is an immediate 
consequence of (121ft . Recall that B is a Brownian motion that starts from x 
under the probability P x , and that I = (I t ,t > 0) is the associated minimum 
process. 

Lemma 2.9 If g is a nonnegative measurable function on R 3 and x > 0, 



dtg(t,I t ,B t fj = 2 dy dz dtq x+z - 2y (t) g(jt,y, z) (24) 



fTb x px roc poo 

/ dy 

'0 Jy 

In particular, if h is a nonnegative measurable function on R 2 , 

pXb x f x /*°° 



E x (J o ° dth{I t ,B t ))=2 jf dy jf dzhfaz). (25) 

Proof of Proposition 12.81 This is a simple consequence of Lemma 12.91 For 
p = 1, the result is exactly formula (|23[1 in Proposition 12.61 We proceed by 
induction on p using the Markov property under n (property (ii) in Proposition 
|2~B|> and then ffifr 

n( / dt 1 ...dt p f(m(t l7 t 2 ) 7 ...,m(t p - l7 t p ),e(t 1 ) 7 ...,e(t p ))] 

V i{0<ti<-<t p <<r} 7 

(i^i . . . dt p — i 

{0<ti<— <t p _i<er} 
To 

dt f(m(t 1 ,t 2 ), ■ ■ . ,m(ip_ 2 ,*p-i),it,e(ii), . . . ,e(t p _i),B t ) J 



e(tp-l) 







2 nl / cfti . . . dtp. 

W{0<tl<-<tp_l<<7} 

dckp_i / <2/3 p f(m{ti,t 2 ), .. , m(i p _ 2 , t p -i), a P -i, e(*i), .. , e(* p _i), /3 P )) . 



The proof is then completed by using the induction hypothesis. □ 

Proof of Theorem 12.71 Let T p be the measurable function from R 2 ^ -1 into 
0( p ) such that 

0(e;*i, . ..,t p ) = Tp(m(ti,t 2 ),.. . , m(t p _i, i p ), e(*i), . . . ,e(t p )). 

The existence of this function easily follows from our construction by induction 
of the marked tree #(e; ti, . . . ,t p ). 
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Denote by A p the measure on R^ -1 denned by 

p-i 

A p (dai . . . dap-xdPi . . . d/3 p ) = (jj l[o,0iA>0 i+1 ]( a i)jdai . . . da p -id/3i . . . d(3 p . 

i=i 

In view of Proposition 12.81 the proof of Theorem 12.71 reduces to checking that 
r p (A p ) = A p . For p = 1, this is obvious. 

Let p > 2 and suppose that the result holds up to order p — 1. For every 
j E {1, . . . ,p — 1}, let Hj be the subset of M^ -1 defined by 

H-j = {(ax, . . . . . . ,P P ); a t > a 3 for every i ^ j}. 

Then, 

p-i 

a p = E 1h i ■ v 

3 = 1 

On the other hand, it is immediate to verify that ■ A p is the image of the 
measure 

Aj(da[ . . . d/3'j) <g> l {0iOo) (h)dh ® A p -j(da" . . . dj3 p _^) 
under the mapping <f> : (a[, . . . ,/3'j, h, a'/ . . . , Pp-j) — ► (oti, ■ ■ • , Pp) defined by 

ctj = h, 

oti = a\ + h for 1 < i < j — 1, 

p i = P i + h for 1 < i < j, 

a, = a"_j +h for j + 1 < i < p— 1, 

A = + ^ for j + 1 < i < p. 

The construction by induction of the tree 9(e; t±, . . . , t p ) exactly shows that, 
a.e. for the measure Aj (da[ . . . dfi'j) ® l^ tOQ )(h)dh ® A p ^j(da'{ . . . d(3p_j), we 
have 

r p o $(a' x , . . . , h, a'{ . . . , = T j (a[, ...,%)* T p ^(a'{ . . . , p'^). 

where if 9 E 0(j) and 9' E ©(p-j), the tree 9 * 9' is obtained by concatenating 

the discrete skeletons of 9 and 9' (as above in the construction by induction of 
9(g; t\, . . . , t p )) and assigning the mark h to the root 0. 

Together with the induction hypothesis, the previous observations imply that 
for any nonnegative measurable function / on 0( p ), 



A p (du)l Hj (u)f(T p (u)) = I dhJJ A,(^')A p -,(du")/(r p ($(u'^,u"))) 

dhff A^du^Ap^du'^f^^^Yp^u")) 



dh / Aj * A p -j(d9)f(6) 

h 
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where we write Aj * A p _j for the image of Aj(dd)A p —j(dd') under the mapping 

h 

(9, 9') — > 9 * 9' . To complete the proof, simply note that 

h 

P-1 poo 
3=1 J ° 

□ 

Remark. The fact that we get only binary trees in Theorem 12.71 corresponds 
to the property that local minima of Brownian motion are distinct: If < t\ < 
■ ■ ■ < t p and if the local minima m g (ti, ti+i) are distinct, the tree 9(g; tt, . . . , t p ) 
is clearly binary. 

2.6. Finite- dimensional marginals of the CRT 

In this subsection, we propose to calculate the law of the tree 9(e; t\, . . . ,t p ) when 
e is a normalized Brownian excursion. This corresponds to choosing p vertices 
independently uniformly on the CRT (the uniform measure on the CRT T e is 
by definition the image of Lebesgue measure on [0, 1] under the mapping p e ) 
and determining the law of the tree spanned by these vertices. In contrast with 
the measure A p of Theorem 12. 71 we will get for every p a probability measure on 
®(pV wn i° n ma y ^ e called the p-dimensional marginal distribution of the CRT 
(cf Aldous [I], 0). 

We first recall the connection between the Ito measure and the normalized 
Brownian excursion. Informally, the law of the normalized Brownian excursion 
(e in the notation of Section^) is n(de | a(e) = 1). More precisely, using a stan- 
dard desintegration theorem for measures, together with the Brownian scaling 
property, one easily shows that there exists a unique collection of probability 
measures (fi( s ), s > 0) on the set £ of excursions, such that the following prop- 
erties hold: 

(i) For every s > 0, n^(a = s) = 1. 

(ii) For every A > and s > 0, the law under n^(de) of e\(t) — ^/\e(t/\) is 

"(As)- 

(iii) For every Borel subset A of £ , 

1 f 00 
n(A) = -(27T)- 1 / 2 J S ^ 2 n (s) (A)d S . 

The measure rt(i) is the law of the normalized Brownian excursion e which 
was considered in Section ^ an d in subsection 12.31 above. 

Our first goal is to get a statement more precise than Theorem 12. 71 by consid- 
ering the pair (9(e; t\, . . . , t p ), a) instead of 9(e; t\, . . . , t p ). If 9 = (r, {h v , v 6 r}) 
is a marked tree, the length of 9 is defined in the obvious way by 
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Proposition 2.10 The law of the pair (9(e; ti,..., t p ), a) under the measure 
n{de) l{o<t 1 < - <t p <a(e)}dti ...dt p 

is 

2 p - 1 A p (d8)q 2L(e) ( S )ds. 



Proof. Recall the notation of the proof of Theorem 12. 71 We will verify that, for 

if, 



any nonnegative measurable function / on 10 '" 



n[ / dti . . . dtp 

' {Q<ti<-<t p <a} 



f(m(ti,t 2 ), ■ ■ ■ ,m(tp-i,tp),e(ti), . . . ,e(t p ),t 1 ,t 2 -ti,...,a- t p )J 
= 2 P ~ 1 / A p (da 1 ..da p - 1 d/3 1 ..df3 p ) / ds x .. ds p+ i q j3l {sx)q^ + p 2 -2 ai (a 2 ) .. 



p+i 



■■ 9/3 p _ 1 +/3 p -2q p _ 1 (sp)q/3 p {s P +i) f{ot\, .. , a p _i,/?i, ...,(3 P , si, .. , s p+1 ). (26) 
Suppose that (|26(l holds. It is easy to check (for instance by induction on p) that 

p-i 

ft, ...,&))= ft + Yj&i + 0i-i - 2a0 + /3 P . 
i=i 

Using the convolution identity q x * q v — q x + y (which is immediate from the 
interpretation of q x as the density of T x under Pq), we get from (|26|l that 



( / dtx... dt p f(m{h, h), m(*p_i, ip), e{h), e(t p ), a) ) 

V J{0<ti< - <t p <a} 1 

2 P ^ 1 A p (d«i . ..dop-idfa . ..d(3 p ) / dt g2i(r„(ai,...,/3 p ))(*) •••,&,*)• 



As in the proof of Theorem 12. 71 the statement of Proposition 12 . 1 01 follows from 
this last identity and the equality r p (A p ) = A p . 

It remains to prove l|26|l ■ The case p = 1 is easy: By using the Markov property 
under the Ito measure fProDOsition l2.6l (ii)). then the definition of the function 
q x and finally (|2"2"|l . we get 

J n{de) J dtf(e(t),t,a-t) = J n(de) J dt E e{t) [f(e(t),t,T )] 



JO 

OO p OO /'OO 



(de) / dt / dt'q e{t) {t')f(e(t),t,t') 



dx / tftg s (t) / dif q x (t') f{x,t,tf). 
o Jo Jo 
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Let p > 2. Applying the Markov property under n successively at t p and at 
tp—i, and then using 1241) . we obtain 



/ dt\ . . . dt p 

J {0<ti<—<t p <a} 

xf(m(ti,t 2 ), ■ ■ .,m(tp- 1 ,t p ),e(t 1 ), . . . , e(t p ), t±, t 2 -ti,...,a- t p )^j 
= n( f dti . . . dtp-i J5 e(fp _ I ) / dt ( dsq Bt (s) 

V J{0<t 1 <---<t p - 1 <cr} l J0 JO 

xf(m(t 1 ,t 2 ), ■■ ,m(tp- 2 ,tp-x),I t ,e(ti), ..,e(tp-i),B t ,ti, .. ,ip_i - t p - 2 ,t,s) J 

= 2n( / dti . . . dtp-i 

^i{0<ii<-<t p _i<o-} 
/■e(i p _i) />oo />oo />oo 

/ dy dz dt ds9 e ( tp _ 1 ) +z _2 3/ (t)g z (s) 

JO Jy JO JO 

x/(m(ti,t 2 ), .. ,m(tp_2,tp-i),2/,e(ti), .. ,e(tp_i), z,ti, .. ,t p _i - tp_ 2 ,t, s)^. 

It is then straightforward to complete the proof by induction on p. □ 

We can now state and prove the main result of this subsection. 
Theorem 2.11 The law of the tree 0{e; ti, . . . , t p ) under the probability measure 

P ] - l{o<ti<-<t„<i}d*i • ■■dtpn(x){de) 



IS 



p\ 2 p+1 L{9) exp ( - 2 L{6) 2 ) A p (dO). 



Proof. Let F be a nonnegative bounded continuous function on ®(p) and let h 
be bounded, nonnegative and measurable on R + . By Proposition ^. 101 



n{de)h{u) I dt x . . . dt p F(9(e;t 1 , . . . ,t p )) 

J{0<t 1 <---<t p <cr} 



= 2"- L J ds h[s) J A p (d6) q 2L{e) {s) F{6). 

On the other hand, using the properties of the definition of the measures nt s \, 
we have also 



n(de)h(a) / dti . . . dt p F(0{e;h, . . . ,t p )) 

J{0<t!<---<t p <(T} 

= i(2 7 r)- 1 / 2 f°° d SS - 3 ' 2 h{s) f n w (de) 



/ 



dt 1 ...dt p F(9(e]t 1 ,...,tp)). 

1 {0<t 1 <-<t p <s} 
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By comparing with the previous identity, we get for a. a. s > 0, 



n( s ){de) I dti...dtpF(8(e;t 1 ,...,tp)) 

J{0<ti<— <t p <s} 

= 2 P+1 J A p (d9) L{9) exp ( - ^^-) F(0). 

Both sides of the previous equality are continuous functions of s (use the scaling 
property of ni s ) for the left side). Thus the equality holds for every s > 0, and 
in particular for s = 1. This completes the proof. □ 

Concluding remarks. If we pick ti, . . . , t p in [0, 1] we can consider the increas- 
ing rearrangement t' x < t' 2 < ■ ■ ■ < t' p of tx, ■ ■ ■ , tp and define 9(e; t\,...,t p ) = 
9(e; t? x , . . . ,t' p ). We can also keep track of the initial ordering and consider the 
tree 9(e; t%, . . . , t p ) defined as the tree 9(e; t\, . . . , t p ) where leaves are labelled 
l,...,p, the leaf corresponding to time ti receiving the label i. (This labelling 
has nothing to do with the ordering of the tree.) Theorem 12 . 1 II implies that the 
law of the tree 9(e; ti, . . . ,t p ) under the probability measure 

l[o,i] p ( < i' ■ ■ ■ >tp)dti ■ ■ .dtpjid^de) 

has density 

2 p+1 L(9) exp(-2L((9) 2 ) 

with respect to Q^{d9), the uniform measure on the set of labelled marked 
(ordered) binary trees with p leaves. 

We can then "forget" the ordering. Let 9(e; t±, .. , t p ) be the tree 9(e; t\, .. , t p ) 
without the order structure. Since there are 2 P ~ 1 possible orderings for a given 
labelled binary tree with p leaves, we get that the law (under the same measure) 
of the tree 9(e; ti,...,t p ) has density 

2 2p L{9) exp(-2L(<9) 2 ) 

with respect to <d^(d9), the uniform measure on the set of labelled marked 
(unordered) binary trees with p leaves. 

In agreement with Aldous' normalization of the CRT, replace the excursion 
e by 2e (this simply means that all marks h v are multiplied by 2). We obtain 
that the law of the tree 9(2e; tx,...,t p ) has density 

with respect to Q^(d9). It is remarkable that the previous density (apparently) 
does not depend on p. 

In the previous form, we recognize the finite-dimensional marginals of Aldous' 
continuum random tree , [5] ■ To give a more explicit description, the discrete 
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skeleton r(2e;ti, . . . ,t p ) is distributed uniformly on the set of labelled rooted 
binary trees with p leaves. This set has b p elements, with 



Then, conditionally on the discrete skeleton, the marks h v are distributed with 
the density 



Bibliographical notes. The coding of real trees described in subsection 12.21 is 
taken from although the underlying ideas can be found in earlier papers 
(see in particular [2] and [53]). A simple approach to Theorem 12.21 based on 
the four-point condition for real trees, can be found in Lemma 3.1 of ^HJ. See 
e.g. Chapter XII of |2U or the last section of for a thorough discussion of 
the Ito excursion measure. The CRT was introduced and studied by Aldous pQ , 
|21 ■ The direct approach to finite-dimensional marginals of the CRT which is 
presented in subsections 12 . 51 and 12 .61 above is taken from |26| . 

3. The Brownian Snake and its Connections 
with Partial Differential Equations 

Our goal in this section is to combine the continuous tree structure studied 
in the previous section with independent spatial motions: In additional to the 
genealogical structure, "individuals" move in space according to the law of a cer- 
tain Markov process £. This motivates the definition of the path- valued process 
called the Brownian snake. We study basic properties of the Brownian snake, 
and we use our previous calculation of marginal distributions of random real 
trees ( subsection 12. 5JI to give explicit formulas for moment functionals. We then 
introduce the exit measure of the Brownian snake from a domain, and we derive 
a key integral equation for the Laplace functional of this random measure. In the 
case when the spatial motion £ is Brownian motion in M. d , this integral equation 
leads to important connections with scmilincar partial differential equations, 
which have been studied recently by several authors. 

3.1. Combining the branching structure of a real tree 
with a spatial displacement 

We consider a Markov process ((,t,Iix)t>o,x£E with values in a Polish space E. 
We will assume that £ has continuous sample paths and that there exist an 
integer m > 2 and positive constants C and e such that for every x £ E and 
t > 0, 



P 



[2-^-^Cp = 1 x 3 x ••■ x (2p-3). 




) 




0<r<t 



(27) 
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where <5 denotes the distance on E. This assumption is not really necessary, but 
it will simplify our treatment. It holds for Brownian motion or for solutions of 
stochastic differential equations with smooth coefficients in R d or on a manifold. 
Later £ will simply be d-dimensional Brownian motion, but for the moment it 
is preferable to argue in a more general setting. 

We denote by W the set of all finite -E-valued paths. An clement of W is a 
continuous mapping w : [0, (] — > E, where £ = C(w) — depends on w and is 
called the lifetime of w. The final point of w will be denoted by w = w((). If 
x G E, the trivial path w such that Q w ) = and w(0) = x is identified with 
the point x, so that E is embedded in W. The set W is a Polish space for the 
distance 

d(w, w') = |C - C'l + sup 5(w(t A C), w'(t A (')). 
t>o 

For i6i?we denote by W x the set W x = {w e W : w(0) = x}. 

If g : M4. — > R+ is a continuous function with compact support such that 
g(0) = 0, we saw in subsection 12 . 21 that g codes a real tree T g . Our goal is now 
to combine this branching structure with spatial motions distributed according 
to the law of the process £. To this end, we will not make an explicit use of the 
tree T g but rather give our definitions in terms of the coding function. It will 
be convenient to drop the compact support assumption on g, and to consider 
instead a general function / 6 C(R+ , R+ ) . 

Notation. Let w : [0, £] — ► E be an element of W, let a £ [0, Q and b> a. We 
denote by i? a) b(w, dw') the unique probability measure on W such that 

(i) C = b, i? a ,fc(w, dw') a.s. 

(ii) w'(t) = w(<), Vi < a, R a _b{w,dw') a.s. 

(iii) The law under i? a ,h(w, dw') of (w'(a + t),0 < t < b — a) coincides with the 
law of (£t , < t < b — a) under II w ( a ) . 

Under i? Qi h(w, dw'), the path w' is the same as w up to time a and then behaves 
according to the spatial motion £ up to time b. 

Let (W s ,s > 0) denote the canonical process on the space C(IR+,W) of 
continuous functions from M. + into W. We also denote by £ s = C(w s ) the lifetime 
ofW s . 

Proposition 3.1 Let f £ C(K+,M+). Assume that f is locally Holder contin- 
uous with exponent r\ for every r] € (0, ^). Let w G W with Q w ) = /(0). Then, 
there exists a unique probability measure on C(R+, W) such that Wq = w ; F{, 
a.s., and, under T^, the canonical process (W s ,s > 0) is (time-inhomogeneous) 
Markov with transition kernel between times s and s' given by 

R mf {s,s')j(s>){W s ,d-w') 

where mf(s, s') = inf r s _ s /j f(r). We have in particular £ s = /(s) for every s > 0, 
Tl a.s. 
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The intuitive meaning of this construction should be clear : At least when 
/(0) = and / has compact support (so that we can use the theory of subsec- 
tion l2.2f> . the vertices p/(s) and Pf(s') in the tree 7} have the same ancestors up 
to generation m/(s,s'). Therefore, the corresponding spatial motions W s and 
W s i must be the same up to time m/(s,s') and then behave independently. 
This means that the path W s < is obtained from the path W s through the kernel 

Rm f (s,s>)J(s')(W s , dw'). 

Proof. We give the detailed argument only in the case when /(0) = 0, which 
implies that w = x G E. We leave the general case as an exercise for the reader 
(cf the proof of Proposition IV. 5 in |27|). 

For each choice of < t\ < t 2 < • • ■ < t p , we can consider the probability 
measure ^tf,...,t p 011 W p defined by 

7T u f .-,t P ( dwi --- dw p'> 

= R Q j( tl ){x, dwi)i? m/(tljt2)i/( t 2) (wi, dw 2 ) . . . i2 m/ ( tp _ 1 ,t f ,),/(i p )(w p _i, dw p ). 

It is easy to verify that this collection is consistent when p and t%, . . . , t p vary. 
Hence the Kolmogorov extension theorem yields the existence of a process 
(W s , s > 0) with values in W (in fact in W x ) whose finite-dimensional marginals 
are the measures 7r? + . 

We then verify that (W s , s > 0) has a continuous modification. Thanks to the 
classical Kolmogorov lemma, it is enough to show that, for every T > there 
are constants f3 > and C such that 

E[d(W s ,W s ,) m ] < C\s-s'\ 1+ P, (28) 

for every s < s' < T. (Here m is as in assumption 127|l .) 

Our assumption on / guarantees that for every -q £ (0, 1/2) there exists a 
finite constant C Vj t such that, for every s, s' £ [0, T], 

\f(s) - f(s')\ < Cr,, T \s - sT- (29) 

By our construction, the joint distribution of (W s , W s i) is 

Roj(s)(x, dw)R mf ( s>s >)j( s >-) (w, dw'). 

This means that W s and W s > are two random paths that coincide up to time 
mj(s,s') and then behave independently according to the law of the process £. 
Using the definition of the distance d, we get for every s, s' € [0, T], s < s', 

E(d(W s ,W s 'D 

<c m (\f( S )-f(s')\ m + 2IL x (u c ( sup ^o,6) m ))) 

0<t<(J(s)Vf(s'))-m f (s,s') 77 

< c m (|/( S ) - f(s')\ m + 2C|(/(s) V f(s')) - m f (s, s ')| 2+£ ) 

< c m (c™ T \s - s'\ mn + ZCCffls - s '|( 2+£ )") , 
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where we used assumption (|27f) in the second inequality, and (|29(l in the third 
one. We can choose r\ close enough to \ so that mr\ > 1 and (2 + e)r\ > 1. The 
bound l(2*51) then follows. 

We then define r£ as the law on C(R+,W) of the continuous modification 
of (W s ,s > 0). The fact that under T* x the canonical process is Markov with 
the given transition kernels is obvious from the choice of the finite-dimensional 
marginals. The form of the marginals also readily shows that Q tl = f{t\), • • ■ , 
Ct p = f(t p ), r| a.s., for every finite collection of times ti, . . . ,t p . The last asser- 
tion of the proposition then follows from a continuity argument. □ 

The process (W s , s > 0) under the probability measure is sometimes called 
the snake driven by the function / (with spatial motion £ and initial value w). 
From the form of the transition kernels, and an easy continuity argument we 
have, r{, a.s. for every < s < s' , 

W s (t) = W s >{t) for every t < m f (s, s'). 

We sometimes refer to this property as the snake property. Note in particular 
that if x = w(0) we have W s (0) = x for every s > 0, r{, a.s. 



3.2. The Brownian snake 



We now randomize / in the construction of the previous subsection. For every 
r > 0, we denote by P r {df) the law of reflected Brownian motion started at r 
(the law of (\B S \, s > 0) if B is a linear Brownian motion started at r). Then 
P r (df) is a probability measure on the set C(K+,M+). Note that the assumption 
of Proposition ETH holds P r {df) a.s. 

For every s > 0, we denote by p r s (dadb) the law under P r of the pair 

( inf /(«),/(*)). 

0<u<s 

The reflection principle (cf subsection I2.4f> easily gives the explicit form of 
p r s {da, db): 

2(r + b-2a) / (r + b - 2a) 2 ^ 
p s (da,db) = (27rs3)i/2 exp { 2s / ( 0<a<bAr ) 

{r + by- 



+ 2 (2ns)- 1 / 2 exp ( - £± r ) l (Q<b) S (da) db. 



Theorem 3.2 For every w £ W, denote by P w be the probability measure on 
C(R+,W) defined by 

Pw(^)= / P c ^(df)Tl(du;). 



The process (W s , P w )s>o,wew is a continuous (time-homogeneous) Markov pro- 
cess with values in W, with transition kernels 

Q s (wi,dw 2 ) = / Ps (wi) (rfa, d6) Ra,b{wi, dw 2 ). 
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This process is called the Brownian snake with spatial motion £. 

Proof. The continuity is obvious and it is clear that Wo = w, P w a.s. The 
semigroup property for the kernels Q s is also easy to verify. As for the Markov 
property, we write 



E w [F(W sl ,...,W Sp )G(W Sp+1 )} 

= / p C ( w)( d /) / RmAO.uMMfadwi) 

J JWP+ 1 



2(p+l) 

+ 



/( Sp+1 )(w p ,dw p+ i) F(w x , . . . ,w p )G(w p+ i) 
p Q s ^ ) (da u db 1 )p b s \_ Si {da2,db 2 ) . . . p b / p+1 _ ap (da p+ x,db p+ i) 



^ai,6i(w,ciwi) . ..R ap+ubp+1 (wp,dvf p+ i)F(wi, . . . , Wp)G(w p +i) 

w p+1 

( P +l) Psi W '( da l' d& l)Ps2- S i( da 2:^2) ■ ■■p b s P ~- Sp _ 1 (d a P' db p) 

fioi,6i(w,dwi) . . . J R ap!6p (wp_i,dwp)i ;l (wi, . . . ,w p )Q Sp+1 _ Sp G(wp) 
- E W [F(W S1) . . . , W Sp )Q Sp+1 - Sp G(W Sp )]. 
This completes the proof. □ 

Under P w , the process (( s ,s > 0) is a reflected Brownian motion started at 
C(w)- This property is obvious from the last assertion of Proposition IXT1 and the 
very definition of P w - 

The snake property can then be stated in the form 

W a (t) = W s >(t) for all i< inf ( r , 

s<r<s' 

for every s < s', P w a.s. 

In particular, if w(0) = x we have W s € W x for every s > 0, P w a.s. 

We now state the strong Markov property of W, which is very useful in 
applications. We denote by J- s the canonical filtration on C(M.+ ,W) {T s is the 
cr-field generated by W r , < r < s) and as usual we take 



r>s 



Theorem 3.3 The process (W s , P w )s>o,wew *s strong Markov with respect to 
the filtration {J : s +). 

Proof. Let T be a stopping time of the filtration {T s +) such that T < K for 
some K < oo. Let F be bounded and Tt+ measurable, and let $ be a bounded 
measurable function on W. It is enough to prove that for every s > 0, 

E W (F*(W T+S )) = E w (FE Wt (V(W s ))). 
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We may assume that is continuous. Then, 

oo 

E W (F^(W T+S )) = hm 

F*(Wfc±i +5 )) 

oo 

FE Wh±L (*(W s ))). 



n — >oo ' 

fc=0 



In the first equality, we used the continuity of paths and in the second one the 
ordinary Markov property, together with the fact that l{k/ n <T<(k+i)/n} F 1S 
^(fc+ij/n-measurable. At this point, we need an extra argument. We claim that 

limf sup |Ew r (*(W s )) -E Wt (V{W s ))\) =0, P w a.s. (30) 

e l° v t<K,t<r<t+e 1 

Clearly, the desired result follows from (|3(J[1 . because on the set {k/n < T < 
(k + l)/n} we can bound 

\M W (*(Ws))-E Wt (9(W s ))\ < sup \E Wr (y(W s ))-E Wt (V(W s ))\. 

t<K,t<r<t+± 

To prove H3()|) . we write down explicitly 

E Wr (V(W s )) = J P i r (da,db) J R a , b (W r ,dw')V(w'), 

and a similar expression holds for Ew t (^>(W S )) . Set 

c(e) = sup |Cr-Ct| 

t<K,t<r<t+e 

and note that c(e) tends to as e — > 0, P w a.s. Then observe that if t < K and 
t < r < r + e, the paths W r and Wt coincide at least up to time — c ( £ ))+- 
Therefore we have 

R a<b (W r ,dw') = R a , b (W t ,dw') 

for every a < (Ct — c { £ ))+ an( l b > a. The claim H30[) follows from this observation 
and the known explicit form of p^ r (da, db). □ 

Remark. The strong Markov property holds for W even if the underlying spa- 
tial motion £ is not strong Markov. 

3.3. Excursion measures of the Brownian snake 

For every x € E, the excursion measure N x is the cr-finite measure on C(K+, W) 
defined by 

N x (duj) = [ n(de) rS(dw), (31) 
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where n(de) denotes Ito's excursion measure as in Section In particular, 
(Csi s > 0) is distributed under N K according to the Ito measure n(de). As 
in Section [21 we will use the notation a for the duration of the excursion under 
N x (duj): 

a = inf{s > : C s = 0}. 

Note that W s £ W x for every s > 0, N x a.e. 

Let us comment on the intuitive interpretation of N x . As we saw in Section|21 
the excursion (£ s ,0 < s < a) codes a random real tree. If a is a vertex of this 
real tree, corresponding say to s £ [0, a], the path W s gives the spatial positions 
along the line of ancestors of a in the tree, and in particular its terminal point 
W s is the position of a. The snake property reflects the fact that the ancestors of 
the vertices corresponding to s and s' are the same up to generation inf^ s /j £ r . 
To summarize, the paths 

W s , < s < a 

form under N x a "tree of paths" of the spatial motion £, whose genealogical 
structure is governed by the tree coded by the Brownian excursion (£ s , < s < 
a). 

Remark. We know that the process (W S) P w ) s >o, w ew , where the driving ran- 
dom function is reflected Brownian motion, is a continuous strong Markov pro- 
cess. Furthermore, every point x £ E (viewed as an element of W) is regular for 
W, in the sense that P x (T {x} = 0) = 1 if T {x} = inf{.s > 0, W s = x}. This last 
property is immediate from the analogous property for reflected linear Brownian 
motion. Thus it makes sense to consider the excursion measure of W away from 
x, in the sense of the general Ito excursion theory (see e.g. Blumenthal 0]), and 
this excursion measure is easily identified with N x . 

In what follows, since we use the notation (W s ,s > 0) for the canonical 
process on C(R+, W), it is important to realize that the driving random function 
of our Brownian snake (W s , s > 0) is either reflected Brownian motion (under 
the probability measures P w ) or a Brownian excursion (under the excursion 
measures N x ). 

We will make use of the strong Markov property under N K . To state it, it will 
be convenient to introduce the following notation. For every w £ W, we denote 
by the distribution under P w of the process 

w; = w sAa , s > o 

where a := inf{s > : ( s = 0} as above. 

Theorem 3.4 Let T be a stopping time of the filtration (J- s+ ). Assume that 
< T < a, N x a.e. Then, if F and G are nonnegative measurable functionals 
on C(WL+,W), and if F is Tt + -measurable, we have 



N x (f G(W t+s , s > 0)) = m x (f e* Wt (go) 
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If we interpret N x as an excursion measure away from x (as explained in the 
remark above), the preceding theorem becomes a well-known fact of the theory 
of Markov processes: See e.g. 0]. Alternatively, it is also easy to give a direct 
proof of Theorem 13.41 using the same method as in the proof of Theorem 13.31 
together with the simple Markov property under the ltd excursion measure (cf 
Proposition ^. 61 (ii)). 

We will now use Theorem 12.71 to derive some important formulas under the 
excursion measure N x . Let p > 1 be an integer. Recall from subsection 12 . 51 the 
notation 0fc? for the set of all marked trees with p leaves, and let 9 S For 
every x £ E, we associate with 9 a probability measure on W p , denoted by 11^, 
which is defined inductively as follows. 

If p = 1, then 9 = ({0}, h) for some h > and we let U e x = Il£ be the law of 
(6,0 <t <h) under IL^. 

If p > 2, then we can write in a unique way 

9 = 9' * 9" , 

h 

where 9' £ 6j?j", 9" £ Q\"L oV and j £ {1, . . . ,p - 1} (the notation 9' * 6", for 

the concatenation of 8' and 8" with root mark h, was introduced in the proof 
of Theorem II .15|) . We then define Ii e x by 

J H e x (dwi,...,dw p )F(-wi,...,w p ) 

= n*(JJ K( d <> ■ ■ ■ ,^)<(dwi', . . .,dw;_ d ) 

F(£[ ,h] w' 1; . . . , ^ [0 ,fe] vrj,€[o,h] w", . . . , ^ ,h] w p-j)) 

where C[o,ft] w denotes the concatenation (defined in an obvious way) of the 
paths (6,0 <t<h) and (w(t),0< i < C(w))- 

Informally, 11^ is obtained by running independent copies of £ along the 
branches of the tree 9. 

Finally, we recall the notation 6(f; t\, . . . , t p ) from subsection 12. 51 and we let 
A p be as in Theorem 12 . 71 the uniform measure on 

Proposition 3.5 (i) Let f £ C(R+,M+) swc/i that /(0) = 0, and Zei < ti < 
t 2 <---<t p . Then the law under T(. of (W tl , ... ,W tp ) is lf x if ' tu -' tp) . 

(ii) For any nonnegative Borel measurable function F on W p , 



/ ds 1 ...ds p F(W Sl ,... > W 8p ))=2*>- 1 [ A p (de)U e x (F) 

J{0<si<— <s p <<r} ' J 
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Proof. Assertion (i) follows easily from the definition of T{. and the construction 
of the trees 9(f; t\, . . . , t p ). A precise argument can be given using induction on 
p, but we leave details to the reader. To get (ii), we write 



NJ dt 1 ...dt p F(W tl ,...,W ti 

V{o<ti<- •<t P <<j} 

= fn(de) f dt 1 ...dt p T e x (F(W tll ... 1 W tp )) 

J J{Q<t 1 <-<t p <a} V ' 

n{de) [ dt 1 ...dt p Ii e J e;tl >~> t »\F) 

J{0<ti<— <t p <cr} 

= 2"- 1 / A p (de)ILl(F). 

The first equality is the definition of N x , the second one is part (i) of the propo- 
sition, and the last one is Theorem 12. 71 □ 

The cases p = 1 and p — 2 of Proposition ^. 51 (ii) are of special interest. Let us 
rewrite the corresponding formulas in a particular case. Recall that we denote 
by w the terminal point of w. For any nonnegative Borel measurable function g 
on E, we have 



N x (J q dsg(W s )j =n*(J o dtgfo)), 

and 

J dsg{W s )) 2 ) =4n,(| dt ( U ^{J dr ^ r 

Remark. In addition to l|31(l . we could also consider the associated normalized 
excursion measure 

NW(dw)= / n {1) (de)T e x (ckj), (32) 

JC(R + ,R + ) 

where nn) is as in Section[21the law of the normalized Brownian excursion. Let 
Z be the random probability measure on E defined under Nx by 



(Z,g) = f dsg(W s 
Jo 



In the case when £ is Brownian motion in E = R d and x — 0, the random 
measure Z is called ISE (Aldous [3]) for Integrated Super-Brownian Excursion. 
ISE and its variants play an important role in various asymptotics for statistical 
mechanics models (see e.g. [2], |ZI])- In such applications, the explicit formula for 
the moments of the random measure Z is often useful. This formula is proved by 
the same method as Proposition ^ . 51 (ii) . using Theorem l2 . 1 ll instead of Theorem 
12.71 Precisely, we have the following result. 
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Proposition 3.6 For any nonnegative Borel measurable function g on E, 

^ ] {{z,g) p ) 

= P \2 p - 1 J A p {d0)L(9)cxp{-2L(6) 2 ) J Il£(dwi . . . dw p ) gfa) • • • g{w p ). 

We conclude this subsection with an important technical lemma. We fix w G 
W x with Q-w) > and now consider the Brownian snake under P w , that is 
the driving random function of the snake is reflected Brownian motion (and no 
longer a Brownian excursion as in the beginning of this section). We again use 
the notation 

o~ = w£{s > : C = 0} 

and denote by [pn, Pi), i € I the excursion intervals of £ s — inf[o )S i Cr before time 
a. In other words, [pn, /?,), i £ I are the connected components of the open set 
[0, a] n {s > 0, Cs > inf [0 , s ] (A- Then, for every i G I we define W l G C(R+, W) 
by setting for every s > 0, 

= W^+^iU +t), < t < Q ■= C(a I+S )A/3, - U ■ 

From the snake property we have in fact W % G C(R+, W w r^ a .)). 
Lemma 3.7 The "point measure 

iei 

is under P w a Poisson point measure on M+ x C(R+, W) wzi/i intensity 

2 l [0>C(w)] (t)dtN w(t) (cL;) . 

Proof. A well-known theorem of Levy (already used in Section states that, 
if (Pt,t> 0) is a linear Brownian motion started at a, the process (3 t — inf [o,t] /3 r 
is a reflected Brownian motion whose local time at is t — > 2(a — inf[o,t] A-)- 
From this and excursion theory, it follows that the point measure 

iG/ 

is under P w a Poisson point measure with intensity 

2 l [0)Cw ](t)dtra(de) . 

It remains to combine this result with the spatial displacements. 

To this end, fix a function / G C(R + ,M+) such that /(0) = C(w), v{f) = 
inf{i > : f(t) — 0} < oo and / is locally Holder with exponent | — 7 for every 
7 > 0. Recall the notation r£ from subsection 13. II above. Denote by ej,j G J 
the excursions of f(s) — infr 0)S ] f(r) away from before time cr(/), by (a,j,bj), 
j G J the corresponding time intervals, and define for every J G J 



= W (oj+s)A6 . (/(a,) + 1) , < t < }{{a 3 + s) A 6 3 ) - f(aj) , 
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From the definition of 1^, it is easily verified that the processes W J , j € J are 
independent under r{,, with respective distributions r^^ o ^. 

Let F be a bounded nonnegative measurable function on R + x C (R+ , W) , 
such that F(t,cu) = if sup ( s (uj) < 7, for some 7 > 0. Recall the notation 
P r (df) for the law of reflected Brownian motion started at r. By using the last 
observation and then the beginning of the proof, we get 



(exp-^i^,^)) = f P C(w) (rf/)r4(exp-^F(/(a j ),^)) 
iei jeJ 



cxp-2 j dt j n(de)r e w{t) (l - e- F ^>) 



= cxp-2 f <W> dtN w(t) (l-e- F ^)) . 
Jo 



The third equality is the exponential formula for Poisson measures, and the 
last one is the definition of N x . This completes the proof. □ 



3-4- The exit measure 



Let D be an open set in E and fix x € D. For every w e W x set 

r(w)=inf{te[0,C(w)],w(t)gD} , 

where inf = +00. Define 

£ D - {W s (r(Wy);s > 0,t(W s ) < ex.} , 

so that £ D is the set of all exit points from D of the paths W s , for those paths 
that do exit D. Our goal is to construct N x a.e. a random measure that is in 
some sense uniformly spread over £ D . To avoid trivial cases, we first assume 
that 

n x (3t > 0,6 #D) > . (33) 

We start by constructing a continuous increasing process that increases only on 
the set {s > : t(W s ) = ( s }. 

Proposition 3.8 The formula 

dr l{T(w r )«7<T(n<v)+£} 
defines a continuous increasing process (L®,s > 0), N x a.e. or P w a.s. for any 

w e W x . 



L? =hm- / 
e / n 
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Proof. Since N x can be viewed as the excursion measure of W away from x, it 
is enough to prove that the given statement holds under P w . Indeed, if follows 
from the construction of the Ito measure that, for every h > 0, N x (-| sup £ s > h) 
is the law under P^ of the first excursion of W away from x with "height" greater 
than h, and so the result under N x can easily be derived from the case of P x . 

We use the following lemma, where w £ W x is fixed. 

Lemma 3.9 Set "f s = (Cs — t(W s )) + and a s = inf{w > : f£ dr l{ 7r>0 } > s}. 
Then o~ s < oo for every s > 0, P w a.s., and the process T s — 7 CTs is under P w a 
reflected Brownian motion started at (£ w — r(w)) + . 

Proposition !^ . 81 easily follows from Lemma l3~9l Denote by (£ S) s > 0) the local 
time at of T. Then, ¥ x a.s. for every s > 0, 

1 r 

£ s = lim — / dr 1) 

s^a e 

D 



L {0<r r <e}- 

Set A s = J Q S dr l{ 7r >o} and Lf 5 = £a s - We get 

1 f As 1 f s 

L ? = lim 7 / dr 1 {o<r r<£ } = lim- / dr l {0 < 7r < £ } 
6 i- i> £ Jo e t" ^ Jo 



»A 3 

40 e J 1 ~ _J do e J 

which is the desired result. □ 



Proof of Lemma 13.91 For every e > 0, introduce the stopping times 

Sf = infjs > : C, > r(W s ) + e} Tf = infjs > Sf : ( s < t(W s )} 

S? n+1 = mf{s > T* : ( s > t(W s ) + e} T* +1 = inf{ S > S^i ■ Cs < r(W s )} 

We first verify that the stopping times S^ and are finite P w a.s. By apply- 
ing the strong Markov property at inf{s > 0, ( s = 0}, it is enough to consider 
the case when w = x. Still another application of the strong Markov property 
shows that it is enough to verify that Sf < oo a.s. To this end, observe that 
Px(Ci > t(Wi) +e) > (by and because, conditionally on Ci, W\ is a path 
of £ with length £i) and apply the strong Markov property at inf {s > 1, £ s = 0}. 

From the snake property and the continuity of s — > £ s , one easily gets that 
the mapping s — ► 7 S is also continuous. It follows that 75^ = eV (C(w) ~ T ( w )) 
and 7s= = e for n > 2. 

We then claim that, for every n > 1, we have 

T*=M{s>S e n :C=T(W S e)} . 
Indeed the snake property implies that for 



S^<r<mS{s>S^:C = T(W Si )}, 
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the paths W r and Ws^ coincide for t < t(Ws^), so that r(W r ) = t(Ws^)- This 
argument also shows that 7 r = £ r — r(Wgt) for < r <T^. 

From the previous observations and the strong Markov property of the Brow- 
nian snake, we see that the processes 

(7( S e +r ) AT e,r > 0), 71=1,2,... 

are independent and distributed according to the law of a linear Brownian mo- 
tion started at e (at e V (C( w ) — T ( w )) f° r n = 1) and stopped when it hits 0. 
Hence, if 

<r E r =mi\s: / 22 1 [S' ,Tt)(u)du > r L 

the process (7 CT e , r > 0) is obtained by pasting together a linear Brownian motion 
started at e V (Q w ) — r ( w )) an d stopped when it hits 0, with a sequence of 
independent copies of the same process started at e. A simple coupling argument 
shows that (7 CT e , r > 0) converges in distribution as e — > to reflected Brownian 
motion started at (C( w ) — T ( w )) + - The lemma follows since it is clear that o~ e r J, ay 
a.s. for every r > 0. □ 

Definition. TTie exit measure Z D from D is defined under N x by the formula 

(Z D ,g) = [° dL?g(W s ), 
Jo 



From Proposition 13.81 it is easy to obtain that increases only on the 
(closed) set {s £ [0,er] : ( s = t(W s )}. It follows that Z D is (N^ a.e.) supported 
on£ D . 

Let us consider the case when (|33f) does not hold. Then a first moment cal- 
culation using the case p — 1 of Proposition 13. 51 (ii) shows that 

/■OO 

/ ds l{ T (w,)<oo} = , a.e. 
Jo 

Therefore the result of Proposition 13.81 still holds under N x with = for 
every s > 0. Consequently, we take Z D = in that case. 

We will need a first moment formula for L D . With a slight abuse of notation, 
we also denote by r the first exit time from D for £. 

Proposition 3.10 Let denote the law o/(£ r ,0 < r < t) under the subprob- 
ability measure H x (- H {r < oo}) (TI^ is viewed as a measure on W x )- Then, 
for every bounded nonnegative measurable function G on W x , 

N x (£ dL?G(W s )) =Hf(G). 
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In particular, for every bounded nonnegative measurable function g on E, 

N x ((Z D ,g)) =n x (l {T<oo} g(£ r )) . 

Proof. We may assume that G is continuous and bounded, and G(w) = if 
C(w) < K^ 1 or Qr w ) > K, for some K > 0. By Proposition |^| 

f dLf G(W S ) = lim- f dsl {T{Ws)<Cs<T{We)+e} G(W s ) (34) 
Jo £ Jo 

N x a.e. If we can justify the fact that the convergence (|34|1 also holds in L 1 (N a; ), 
we will get from the case p = 1 of Proposition ^. 51 (ii): 

NjJ^ dL?G(W s ?) = linii| dhIl x (l {T<h<T+e} G(Z r ,0<r <h) 

= limn^li^ooje- 1 ^ dhG(£ r ,0<r< /i)) 

r <oo} 



It remains to justify the convergence in L^PLj). Because of our assumption on 
G we may deal with the finite measure n {sup£ s > .ff -1 }) and so it is 
enough to prove that 

sup (W 3 )<( 3 <t(W s )+e} G(W S )) 

0<e<l u e Jo ' ' 

is finite. This easily follows from the case p — 2 of Proposition 13.51 (ii), using 
now the fact that G(w) = if C( w ) > K. □ 

Let us give an important remark. Without any additional effort, the previous 
construction applies to the more general case of a space-time open set D C 
K + x E, such that (0, x) G D. In this setting, Z D is a random measure on 
dD cM+x E such that for g G C b+ {dD) 

1 f' 7 

(2 D ,ff) = Hm- / ds l{ T {w s )<t s <T{w s )+e}9(C,s,W s ) 
£ Jo 

where t(w) = infjt > 0, (t,w(t)) ^ D}. To see that this more general case is 
in fact contained in the previous construction, simply replace £ by the space- 
time process = (i,£t); which also satisfies assumption l|2T|). and note that the 
Brownian snake with spatial motion £' is related to the Brownian snake with 
spatial motion £ in a trivial manner. 



We will now derive an integral equation for the Laplace functional of the 
exit measure. This result is the key to the connections with partial differential 
equations that will be investigated later. 
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Theorem 3.11 Let g be a nonnegative bounded measurable function on E. For 
every x € E, set 

u(x) = N x (l- exp -(Z D ,g)) , xeD. 
The function u solves the integral equation 



u{x) + 2JL t (^J u^ s ) 2 ds^ =U x (l {T<oo} g(^)) 



(35) 



Our proof of Theorem 13.111 is based on Lemma 13.71 Another more computa- 
tional proof would rely on calculations of moments of the exit measure from 
Proposition 13. 51 above. 



Proof. For every r > set rff = inf{s > : > r}, with the usual convention 
inf = oo. By the definition of Z D , we have 



u(x) = N x 

= N x 



1 — exp - 



dL?g(W s ) 



dLfg(W s )exp(- / dL?g(W r )) 



drl {v?«*>}9( W v?) exp- / dL s g(W, 



dr 1{ v d <co 



dL°g(W s )E w Jexp- I dL?g(W r ) 



} g(W n o ) Ew ? (exp - f dL?g(W s ) 



The second equality is the simple identity 1 — exp(— A t ) = f Q dA s exp(— (A t — 
A s )) valid for any continuous nondecreasing function A. The third equality is the 
change of variables s — rff and the fourth one follows from the strong Markov 
property under (cf Theorem 13.4(1 at the stopping time rff . 

Let w £ yV x be such that £( w ) = t(w). From Lemma 13.71 we have 



E w (exp-^ dL?g(W r 



E w ( exp — ^ 

/•C(w) 

exp 



Pi 



dL°g(W r ) 



(i 



exp I —2 



dtu(w(t)) 



Hence, 



u{x) = E X (J^ dL?g(W s )exp(-2 J' dtu(W s {t)))) 
= tt x (l {T<oc} g(£ T )exp(-2 J dtufa)) 
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by Proposition y . 1 ()1 The proof is now easily completed by the usual Feynman- 
Kac argument: 

u(x) = U x (l {T<oc} g(^ T )) - n a; ^l {r<oo}i g(^ 7 -)(l - exp-2 J dt 

= n,(l {T<oo} .g(e r )) -2n.(l {T<oo} .g(e r )^ dtu(^)exp(-2^ dr «&.))) 

= n x (1{ t<0 o}5(Ct)) - 2n a ^ dbu(£t)TLt t (l{ T<00 yg(£ T )exp(-2j dru(£ 

= (l {T<0 o}5(er )) - 2H X ( / dt M (6) 2 ) • 

Jo 

3.5. TTie probabilistic solution of the nonlinear Dirichlet problem 

In this subsection, we assume that £ is Brownian motion in M d . The results 
however could easily be extended to an elliptic diffusion process in M. d or on a 
manifold. 

We say that y e dD is regular for D c if 

inf{f > : & g £>} = , IT, a.s. 

The open set D is called regular if every point y £ dD is regular for D c . We say 
that a real- valued function u defined on D solves Am = 4u 2 in D if u is of class 
C 2 on £> and the equality Au = Au 2 holds pointwise on D. 

Theorem 3.12 Let D be a domain in M. d and let g be a bounded nonnegative 
measurable function on dD. For every x € D, setu(x) = N x (l-exp-(Z D ,g)). 
Then u solves Au = 4m 2 in D. If in addition D is regular and g is continuous, 
then u solves the problem 

Au = 4m 2 in D 

(36) 

u\dD =9 

where the notation u\qd = g means that for every y £ dD, 

lim u(x) = g{y) . 

L)3x — >y 

Proof. First observe that, by i|35[l . 

u(x) < IL c (li T<ao \g(£ T j) < sup g(y) , 

yedD 

so that u is bounded in D. Let B be a ball whose closure is contained in D, 
and denote by tb the first exit time from B. From Q35|l and the strong Markov 
property at time tb we get for x G B 



u{x) + 2U X U B u(£ s ) 2 ds) + 2U X (ll iTB ( J 
= U x (U^ B (l {T<oo} g^ r ))). 



u(£, s ) 2 ds 
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By combining this with formula (|35[1 applied with x — £ TB , we arrive at 

u{x) + 2IL X (y B u^ s ) 2 ds^j = Tl x (u^ TB )) . (37) 

The function h(x) = U x (u(£ TB )) is harmonic in B, so that h is of class C 2 and 
Ah = in B. Set 

/(*) :=UJ[ TB u(0 2 ds) = f dyG B (x,y)u(y) 2 

JO J B 

where Gb is the Green function of Brownian motion in B. Since u is measurable 
and bounded, Theorem 6.6 of shows that / is continuously differentiable in 
B, and so is u since u = h—2f. Then again by Theorem 6.6 of |37j . the previous 
formula for / implies that / is of class C 2 in B and —A Af = u 2 in B, which 
leads to the desired equation for u. 

For the second part of the theorem, suppose first that D is bounded, and let 
y S dD be regular for D c . Then, if g is continuous at y, it is well known that 

lim n x (g(£ T )) = g(y) . 
On the other hand, we have also 

limsupn^f / u(t; s ) 2 ds) < (sup u(x)) 2 lim sup E x (r) = . 



Thus l|33|l implies that 



lim u(x) = g(y) 

D3x — >y 



When D is unbounded, a similar argument applies after replacing D by D n B, 
where B is now a large ball: Argue as in the derivation of Ij37(l to verify that for 
x £ DC\B, 

u(x) + 2U x ^ u(^) 2 ds) =U x (l {T < TB} g{^)) +U x (l {TB<T} u(^ TB )) 

and then follow the same route as in the bounded case. □ 
The nonnegative solution of the problem l|36l) is always unique. When D is 
bounded, this is a consequence of the following analytic lemma. 

Lemma 3.13 (Comparison principle) Let h : R+ — > R+ be a monotone increas- 
ing function. Let D be a bounded domain in R d and let u, v be two nonnegative 
functions of class C 2 on D such that Au > h{u) and Av < h(v). Suppose that 
for every y £ dD, 

limsup(w(x) — v(x)) < . 

DBx — >y 

Then u < v. 
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Proof. Set / = u — v and D' = {x £ D, f(x) > 0}. If D' is not empty, we have 
A/0) > /i(u(a;)) - h(v(x)) > 

for every x £ D' . Furthermore, it follows from the assumption and the definition 
of D' that 

limsup f{x) < 

for every z S Then the classical maximum principle implies that / < on 
D' , which is a contradiction. □ 

Corollary 3.14 (Mean value property) Let D be a domain in R d and let U be 

a bounded regular subdomain of D whose closure is contained in D. Then, if u 
is a nonnegative solution of Au = 4 u 2 in D, we have for every x G U 

u{x) =N x (l-exp-(Z c/ ,u)). 

Proof. For every x € U, set 

v{x) =N x (l-cxp-(Z u ,u)). 

By Theorem 13.121 v solves Au = Av 2 in U with boundary value v\qjj — u\qu. 
By Lemma 131 we must have v(x) — u(x) for every x G U . □ 

The last proposition of this subsection provides some useful properties of 
nonnegative solutions of Au = 4u 2 in a domain. For x £ M. d and e > 0, we 
denote by B(x,e) the open ball of radius e centered at x. We also denote by 

K = {W s ,0 < s<a} 

the range of the Brownian snake. 

Proposition 3.15 (i) There exists a positive constant Cd such that for every 
x e K d and e > 0, 

n x {nr\B{ x ,e) c ^0) =c d£ - 2 . 

(ii) Let u be a nonnegative solution of Au = Au 2 in the domain D. Then for 
every x G D. 

u{x) < c d dist(x,dDy 2 . 

(iii) The set of all nonnegative solutions of Au = Au 2 in D is closed under 
pointwise convergence. 

Proof, (i) By translation invariance we may assume that x = 0. We then use a 
scaling argument. For A > 0, the law under n(de) of e\(s) = A _1 e(A 2 s) is A _1 n 
(exercise !). It easily follows that the law under No of Ws 6 \t) = E~ 1 W e i s (e 2 t) 
is £~ 2 No- Then, with an obvious notation, 

n a (nr\B(o,e) c ^ 0) = No(n ie) ns(o, i) c ^ 0) 



= £- 2 N (^ns(o,i) c ^ 0) 
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It remains to verify that No (71 D -8(0, l) c 7^ 0) < 00. If this were not true, 
excursion theory would imply that Po a.s., infinitely many excursions of the 
Brownian snake exit the ball B(0, 1) before time 1. Clearly this would contradict 
the continuity of s — + W s under Po. 

(ii) Let x € D and r > be such that B(x,r) C D. By Corollary 13.141 we 
have for every y £ B(x,r) 

u(y)=N y (l-exp-(Z B ^ r \u)) . 

In particular, 

u(x) < N x (Z B ^ ^0)<K x (TZn B(x, rf ± 0) = c d r~ 2 . 

In the second inequality we used the fact that Z B ^ x,r ^ is supported on E B ( x > r ) c 
1lC\B(x,r) c . 

(iii) Let (u n ) be a sequence of nonnegative solutions of Au = 4u 2 in D such 
that u n (x) — ► u(x) as n — > 00 for every x £ D. Let C/ be an open ball whose 
closure is contained in D. By Corollary 13. 141 for every n > 1 and x E U, 

u n (x) = N x (l - exp - (Z u , tt n )) . 

Note that N a; (-Z' 7 7^ 0) < 00 (by (i)) and that the functions u n are uniformly 
bounded on dU (by (ii)). Hence we can pass to the limit in the previous formula 
and get u(x) — N x (l — exp — (Z u , u)) for x £ U. The desired result then follows 
from Theorem 13.121 □ 

Let us conclude this subsection with the following remark. Theorem 13.111 
could be applied as well to treat parabolic problems for the operator Au — 4m 2 . 
To this end we need only replace the Brownian motion £ by the space-time 
process (i, £t). If we make this replacement and let D C M+ x W 1 be a space- 
time domain, then for every bounded nonnegative measurable function g on dD, 
the formula 

u(t,x) =N t ,,(l-cxp-(Z I? , ff )) 

gives a solution of 

du 1 . „ 9 
— + -Au - 2u 2 = 
9* 2 

in _D. Furthermore, u has boundary condition g under suitable conditions on D 
and g. The proof proceeds from the integral equation l|35|) as for Theorem 13. Ill 

3.6. Solutions with boundary blow-up 

Proposition 3.16 Let D be a bounded regular domain. Thenu\{x) = N X (Z D 7^ 
0), x € D is the minimal nonnegative solution of the problem 

Au = 4u 2 in D 

(38) 
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Proof. First note that u\ (x) < oo by Proposition 13.151 (i) . For every n > 1, 
set v n (x) = N x (l - exp-n{Z D ,l)), x £ D. By Theorem |3~T21 v n solves J2>) 
with g = n. By Proposition 13 . 1 51 (iii), ui = lim j w„ also solves Au = 4m in 
Z?. The condition z*i |az? = 00 is clear since u\> v n and v n \do = n. Finally if 

v is another nonnegative solution of the problem (|38|1 , the comparison principle 
fLemma l3.13ll implies that v > v n for every n and so v > u\. □ 

Proposition 3.17 Let D be any open set in M. d and v,2(x) = N X (7Z n D c ^ 0) 
for x £ D. Then U2 is the maximal nonnegative solution of Au = 4it 2 in D (in 
the sense that u < U2 for any other nonnegative solution u in D). 

Proof First note that 1Z is connected N x a.e. as the range of the continuous 
mapping s — > W s . It follows that we may deal separately with each connected 
component of D, and thus assume that D is a domain. Then we can easily 
construct a sequence (D n ) of bounded regular subdomains of D, such that D = 
lim | D n and D n C D n+ \ for every n. Set 



for x € By the support property of the exit measure, it is clear that v n < v n . 
We also claim that t5 n _|_i(ar) < v n (x) for x € D n . To verify this, observe that on 
the event {Knf^ +1 ^ 0} there exists a path VF S that hits D° +1 . For this path 
W s , we must have TD n (W s ) < ( s (here Tfl n stands for the exit time from D n ), 
and it follows from the properties of the Brownian snake that 



N x a.e. on {TZDD'f l+1 ^ 0}. However, from the construction of the exit measure 
in subsection 13.41 above. (Z Dn ,l) is obtained as the local time at level and 
at time j4" of a reflected Brownian motion started at 0. Since the local time at 
of a reflected Brownian motion started at immediately becomes (strictly) 
positive, it follows that {KflD^ +1 ^ 0} C {Z D ™ ^ 0}, N x a.e., which gives the 
inequality v n +i(x) < v n (x). 

We have then for x € D 



This follows easily from the fact that the event {1Z n D c ^ 0} is equal N x a.e. 
to the intersection of the events {TZD ^ 0}. By Proposition 13 . 1 61 v n solves 
Au = Au 2 in D n . It then follows from H39JI and Proposition 13.151 (iii) that U2 
solves Au = 4u 2 in D. Finally, if u is another nonnegative solution in D, the 
comparison principle implies that u < v n in D n and it follows that u < U2- Q 

Example. Let us apply the previous proposition to compute N K (0 £ TV) for 
x 7^ 0. By rotational invariance and the same scaling argument as in the proof 



v n {x) = N X (Z D " ? 0) , v n (x) = N X (TZ n D c n ± 0) 




u 2 (x) = lim I v n (x) = lim J. v n {x) . 



(39) 
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of Proposition l3~T5l fi). we get N x (0 £ TV) = C\x\~ 2 with a nonnegative constant 
C. On the other hand, by Proposition 13 . 1 7l we know that u(x) = N x (0 € TV) 
solves Aw = 4it 2 in K d \{0}. A short calculation, using the expression of the 
Laplacian for a radial function, shows that the only possible values of C are 
C = and C = 2 — ^. Since u is the maximal solution, we conclude that if 
d < 3, 

N x (0eK) = (2-±)\x\- 2 

whereas N x (0 € TV) = if d > 4. In particular, points are polar (in the sense 
that they are not hit by the range) if and only if d > 4. 

Let us conclude with some remarks. First note that, if D is bounded and regular 
(the boundedness is superfluous here), the function u-i of Proposition 2 also 
satisfies i^lau = +oo. This is obvious since Ui > u%. We may ask the following 
two questions. 

1. If D is regular, is it true that u\ — u-p- (uniqueness of the solution with 
boundary blow-up) 

2. For a general domain D, when is it true that U2\od = +oo? (existence of 
a solution with boundary blow-up) 

A complete answer to question 2 is provided in [H] (see also [231 ) - A general 
answer to 1 is still an open problem (see and the references therein for 
partial results). 

Bibliographical notes. Much of this section is taken from |27j . where additional 
references about the Brownian snake can be found. The connections with partial 
differential equations that are discussed in subsect ions !3 . 5l and l3 . 61 were originally 
formulated by Dynkin \Y6\ in the language of superprocesses (see Perkins |35j 
for a recent account of the theory of superprocesses). These connections are 
still the subject of an active research: See Dynkin's books E3- Mselati's 
thesis [HH] gives an application of the Brownian snake to the classification and 
probabilistic representation of the solutions of Au — u 2 in a smooth domain. 
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